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investigators. In figure 6 the logarithm law is verified. 
Figure 7 show the variation of the longitudinal fluctua- 
ting velocity with radial distance and comparison is 
made with Lawn results. The values obtained have the 
same relative but slightly higher. This fact agrees with 
the variation of fluctuating velocity at the centre line 
of the pipe with Reynolds number as given by coantic 
[4]. Tse results shown on figures 2 to 7 demonstrate 
that the flow emerging from the pipe is in good 
agreement with previous measurements and corresponds 
to fully developed pipe flow. 


Annulus x = .483 Re=406x 1014 Measurements at 5 mm upstream of the exit the 
+ Horizontal Right annulus are presented on figures 8 to 13. Figure 8 
A » Left Re=1D, —D,)U/y ; 
demonstrates the symmetry of the flow by super-impo- 
Y Vertical Upper . 0 
a Lodo sing measured profiles obtained 90º apart around the 
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annulus. The maximum deviation at the maximum velo- 
city, is less than 1,0%. Figue 9 shows profiles for two 
Reynalds numbers and relates them to the measurements 
of Durst [6]. Figure 10 makes a similar comparison with 
the measurements of Brighton and Jones [2]. The 
annular flow has two important parameters to be consi- 
dered: Reynolds number and the inner to outer diameter 
ratio. The comparison became then more difficult but 
in spite of this it can be seen that the location of the 
maximum velocity is consistant with the results presen- 
ted by other investigators, The mean velocity curves in 
general have very close agreement. Figure 11 presents 
profiles for two Reynolds numbers and for the inner and 
outer wall of the annulus on ut — y* co-ordinates and 
comparison with the measurements of Quarmby [15] 
and Brighton and Jones is made in Figures 12, The 
variation of the longitudinal fluctuating velocity with 
radial distance is shown in figure 13 and comparison 
is made with the results of Lawn and Elliott [13] and 
Brighton and Jones. Again, the results indicate that the 
flow is in good agreement with measurements obtained 
by previous workers and that it can safely be considered 
fully developed annulus flow. It is clear from the mea- 
surements that the three pin support does not signifi- 
cantly influence the mean profiles at the annulus exit. 
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3.2 Downstream jet flow fimo 
Yo=1 

The symmetry of the jet is confirmed by figures LO ME, e ! | SJ vertical 
14 to 16 which, for velocity ratios of 0,1.72 and 14,65 pus: ty bo AM Right 
respectively superimpose vertical and horizontal profiles uv” A a ” co Horizontal (rare 
of mean velocity at downstream distances of 10 and 30 rm “Vo % ae * Wertical (a 
diameters. The deviation is less than 2% and may be b Vs o| ' Ei á 
considered very satisfatory. o di id NS 

Figure 17 shows the superposition of mean-velo- 4 by Rays (04) men) 
city measurements at x/D, = 30 (X/D, = 10.86) for x BA 
different Reynolds number ratios and it can be seen that é Ha O Ns 
similarity of mean-velocity profiles has been obtained. AM a 
The figure shows U/U ... as a function of y/x and 4 TENRA a 


allows comparison with the round-jet measurement of 
Hinze [10], Rodi [16], and Wygnanski and Fiedler Fig. 15— U/U max 1Y/Y mo) 
[17]. Measurements in the fully developed region of a 
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Fig. 21 — E and Ve? function of radial distance 


round jet are compared with abramovich [1], champa- 
gne and Wygnanski [3] and Gibson [9] results in fi- 
gure 18. Figure 19 shows the present measurements 
compared with the results of Champagne and Wygnanski; 
the inflence of U, /U, ratio on the spread of the inner 
potential core can be seen; it is also clear that, with a 
higher initial intensity of turbulence, the self preserving 
flow is reached more rapidly than with a lower initial 
intensity of turbulence. (*) 

Hot-wiret transverses made at 30 diameters downs- 
tream for three different Reynolds number ratios are 
reproduced on figures 20 to 22 and again show a very 
satisfatory symmetry. On these figures the signal E 


and e” corresponding to the mean and rms components 
of turbulence indicate a maximum deviation less 
than 2%. 


4 CONCLUSIONS 


(a) Mean velocity measurements at Reynolds num- 
ber from 1.70x104 to 5.08xX104 show that the pipe 
annulus flows are fully developed at the exit. The 
agreement with measurements of other investigators is 
excellent. 


(*) More detailed discussed in reference 18. 
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(b) The maximum discrepancy between measure- 
ments of mean velocity at x/D, equal to O, 10 and 30 
at each of four quadrants is less than 2%. The rms signal 
at x/D, — 30 again measured on each of the four qua- 
drants, also reveals discrepancies less than 2%. The 
concentricity of the two jets is therefore confirmed. 

(c) The mean velocity profiles show that the jet 
is fully developed at x/D, = 30 (x/D, — 10.86). 

(d) The self-preserving jet flow is reached more 
rapidly with higher initial intensities of turbulence. 

(e) The support pins do not significantly influence 
the mean profiles at the annulus exit nor the rms signal 
at a downstream distance equal to 10.86 outer dia- 
meters. 


5. FUTURE WORK 


(a) Hot wire measurements obtained in the free 
jet situation and presented in this report have not been 
evaluated in terms of fluctuating velociies. The turbu- 
lence intensities occuring in the present flow configura- 
tion include values in excess of 20% and necessitate 
approximate means of signal analysis. The means of 
converting the electrical signal into mean fluctuating 
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velocity values is presently being coded according to 
the suggestion of reference 7, 8 and 16. 

(b) The measurements will be extented and detai- 
led investigation of the normal and Reynolds stresses 
will carry on in the co-axial developing region, These 
measurements will assist the development of a turbu- 
lence model. 

Additional measurements with the laser anemome- 
ter will be effected in the present flow configurations. 
These cannot be specified until the above program is 
completed. (*) 
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NOMENCLATURE 

C, = skin friction factor 

D = pipe diameter, 

E=E+e = electrical signal from the hot-wire. 

L = lenght of the pipe. 

R = pipe radius. 

Re = Reynolds number. 

U = mean velocity 

u = fluctuating velocity 

Ur U/U- 

U- V- o = friction velocity 

MD ME - = coordinate system. 

pego % 

Ym2 = distance from the jet axis to the 
radial location where the velocity 
is half the maximum velocity. 

ce R, Ro 


“ 


kinematic viscosity 
ç — density of the flow 
: = wall shear stress 


SUBSCRIPTS 
i = inner side, 
o = outer side. 
m = mean 


max = maximum 


investigation performed in the experimental set up described 


here and pointed out in section 5 has already been, writeten, see references 18, 19 and 20 
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pois irá indicar quando e porque se devem substi- 
tuir. Vai ajudar, também, os nossos concor- 
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DETERMINAÇÃO EXPERIMENTAL DA FUNÇÃO DE TRANSFERÊNCIA 
DO REACTOR PORTUGUÊS DE INVESTIGAÇÃO À POTÊNCIA ZERO 


RESUMO 


Mediu-se a função de transferência do Reactor 
Português de Investigação (RPI), a baixa potência, pelo 
método de modulação da reactividade, na banda de 
frequências compreendida entre 5.10 — 3 Hz e 65 Hz. 
Descrevem-se as características dos dispositivos expe- 
rimentais e da técnica de medição utilizados, bem como 
o princípio do tratamento numérico dos dados experi- 
mentais. 

O modulador concebido para esta experiência pos- 
suía boas características de estabilidade para as frequên- 
cias estudadas e a técnica adoptada na detecção e aná- 
lise da resposta do reactor revelou-se muito eficaz, Os 
resultados obtidos são comparados com os fornecidos 
pela teoria cinética pontual, constatando-se a influência 
dos harmónicos espaciais no comportamento dinâmico 
do reactor, a qual não é tida em conta pelo modelo 
pontual. Dos valores da função de transferência deduz- 
-se para o tempo de geração dos neutrões instantâneos 
no RPl o valor A = (56 + 2).10 6 s para uma frac- 
ção efectiva dos neutrões retardados 5, = (760 + 23). 
ro — 3. 


INTRODUÇÃO 


O comportamento dinâmico de um reactor nuclear 
pode ser descrito pela sua função de transferência 
[1, 2, 3], ou seja, pela resposta a pequenas variações 
sinusoidais de reactividade. A baixa potência, quando o 
combustível e o fluido de arrefecimento não sofrem 
um aquecimento sensível, a função de transferência é 
determinada pelas características da emissão dos neu- 
trões retardados e pelo tempo de geração dos neutrões 
instantâneos. 

O método de modulação de reactividade é um 
processo clássico [1, 3, 7] de determinação experi- 
mental da resposta do reactor (amplitude e fase), fre- 
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ABSTRACT 


The zero power transfer function of the Portuguese 
Research Reactor (RPI) was measured using the reacti- 
vity oscillation technique over the frequency range of 
5.10 —- 3 Hz to 65 Hz. The experimental set-up and 
measurement technique used are described, as well 
as the numerical handling of the data. 

The reactivity oscillator designed for the present 
experiment had good stability at the frequencies studied 
and the technique used in the detection and analysis of 
the reactor output proved to be very efficient. The re- 
sults obtained are compared to the point-reactor kinetics 
model and the effect of the space harmonics on reactor 
dynamics is observed. From the measured transfer func- 
tion a prompt neutron generation time of a = (56 + 2). 
10-6 sec is obtained with an effective delayed neu- 
tron fraction 8, = (760 +23). 10 —S. 


quência por frequência, A excitação periódica é provo- 
cada por um modulador que, no caso presente, foi 
concebido por forma a poder ocupar a posição de um 
elemento de combustível do RPI [4]. 

Neste trabalho são apresentadas as características 
dos dispositivos experimentais e do método de medi- 
ção que foram utilizados na análise da função de trans- 
ferência do RPI a baixa potência, na banda de frequên- 
cias compreendida entre 5.10 —3 Hz e 65 Hz. São 
também indicados o princípio do tratamento numérico 
dos dados experimentais e as correcções que é neces- 
sário introduzir. 

Os resultados obtidos são comparados com a teo- 
ria cinética pontual. 
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1 — PRINCÍPIO 


A análise rigorosa do comportamento dinâmico 
de um reactor exige a consideração da geometria e do 
espectro de neutrões característicos do sistema. É, toda- 
via, usual adoptar-se uma aproximação independente 
do espaço e a um grupo de energias em que se 
admite que: 


— O fluxo de neutrões é uma função exclusiva 
do tempo. 

— À distribuição espacial e energética do fluxo 
não é modificada (coincide sempre com o modo 
fundamental). 


Em geral, este modelo pontual é válido no caso 
de variações de reactividade que se distribuem unifor- 
memente no volume do núcleo ou no caso de variações 
localizadas de reactividade cuja resposta é observada a 
alguma distância (vários livres percursos médios de 
transporte) do local de excitação. 

As equações da cinética tomam então a seguinte 
forma (para não sobrecarregar o texto, o significado dos 
símbolos é indicado no fim): 


dn (t) lt) — B. m 
= n(t+ 2% 4 C tt) (1,1) 
dt A 1-1 
de d(U 8 
em n(t)—> Cat) (1.2) 


Nestas equações faz-se figurar o tempo de geração 4 
(e não o tempo de vida |) como invariante, o que é 
recomendável quando se examinam variações da secção 
eficaz macroscópica de absorção dos neutrões. 

Pretende-se analisar a resposta do reactor, inicial- 
mente crítico ( 2, = O), a oscilações de reactividade 
da fraca amplitude (aproximação linear): 


e()-—- 0, + Aclt) (1.3) 


q 


Seguindo o método clássico de resolução, obtém-se 
a função de transferência do reactor a baixa potên- 
cia [3,6]: 


Nº jo) !n | “ 8 
H (j0) - em joj A+ 4 = 
R (ju) ê > 


H(j vw) é uma função complexa cuja amplitude e fase 
são as da função de transferência do reactor. 

Em particular, se a variação de reactividade for 
sinusoidal 


Ae (t) =: sen (ut), (1.5) 


obtém-se, por inversão de (1.4): 


Antb'n, 
= | H (jo) | sen ot d (1). (1.6) 


E 
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Na prática, a variação de reactividade é normal- 
mente uma soma de harmónicos de período funda- 
mental T: 


op ejhut (1.7) 


BO = E) — 8 = 
| 


h 


[E 


A função de transferência deduzida directamente a 
partir das equações da cinética toma a forma [9]: 


» » x a » na - 1 
np =Hy (fo Mas tea Mp 2 co +tênc Mo) 


com 


E À Zn ===] 
H, — H (ho) jhoj a + 5 E Rm 
jho + 4; 


(1.9) 
Para o fundamental (h = 1) tem-se: 
ny = H, 41h, n, el 1.10) 
ou seja 
a n;/n 
cj =—€ eJri sb? : di eJtyy — Pp), 
H, 
(1.11) 
n,/nm 
H = —— (1.12) 
[é 
"1 
9 
da 4 “1 (1.13) 


Portanto, para obter H, é preciso conhecer <:, ; caso 
contrário, apenas se disporá do valor relativo. Já a fase 
bt, pode ser determinada em valor absoluto. 


2 — EQUIPAMENTO 
2.1 — Modulação da reactividade 


O reactor é excitado quase sinusoidalmente por 
meio de um modulador rotativo com quadrantes de 
cádmio constituído por um rotor cilíndrico que se move 
no interior de um tubo em alumínio (stator). O con- 
junto pode ocupar a posição de um elemento de com- 
bustível na grelha do núcleo do RPI (fig. 1), sendo o 
rotor posto em movimento sob a acção de um motor 
eléctrico, com a potência de 1,5 HP a 1420 rpm, ali- 
mentado com corrente trifásica. O motor e o disposi- 
tivo de transmissão do movimento estão apoiados numa 
plataforma situada na ponte do reactor a 9,8 m da gre- 
lha. Um variador de velocidade acoplado ao motor per- 
mite cobrir uma gama contínua de O a 1420 rpm, em 
ambos os sentidos, 


TÉCNICA 419 


e 


C 
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CC - Camara de cisão 


normal 
Elemento ge combustivel especial 


Cl-camara de tonização 


Fig. 1 — Configuração do núcleo do R.P,l. (14-1) 


Sobre o stator e o rotor, a meia altura do núcleo, 
estavam montados rectângulos de cádmio com dimen- 
sões ajustadas às condições experimentais. Assim, para 
provocar quatro períodos de excitação por período do 
movimento, instalavam-se no rotor quatro rectângulos 
de cádmio com altura h= 25 mm e base b= p/8 
(p = perímetro do rotor); no stator fixavam-se também 
quatro rectângulos com h' = h+ 2 eb'= p'/8 (p' = pe- 
rimetro do stator). Para provocar dois períodos de excita- 
ção, montavam-se dois rectângulos com h = 15 mm e 

= p/4 no rotor e outros dois com h'=h+2 e 
b' =p'/4 no stator. Para provocar um período de excita- 
ção por período do movimento: um rectângulo com 
h=13mmeb=p/2 norotore outro com h'=h+2 
e b'= p'/2 no stator. 

Às gamas de frequência estudadas foram as se- 
guintes: 


— Com 4 períodos de excitação/volta: 3 a 65 Hz. 

— Com 2 períodos de excitação/volta: ta 33 Hz. 

— Com 1 período de excitação/volta: 0,5 a 
16,5 Hz. 
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Os valores indicados dizem respeito às experiên- 
cias efectuadas com o motor ligado directamente ao 
velo do rotor, Para os estudos a frequências mais bai- 
xas, dispunha-sa de um elstema de desmultiplicação 
com «duas velocidades» que permitia cobrir uma gama 
de frequência entre 5.10 — 3 6 2Hz (com um período 
de excitação/volta). 


2.2 — Medição do periodo da modulação 


A medição do período da modulação faz-se atra- 
vés da contagem dos impulsos fornecidos por um gera- 
dor funcionando a uma frequência fixa e conhecida 
(1 kHz). O acesso ao contador é condicionado pelo 
estado de uma porta que permanece aberta durante um 
número inteiro de períodos de rotação previamente fixado. 
O funcionamento desta porta é comandado por impulsos 
que são fornecidos por uma foto-resistência de cada 
vez que esta deixa de ser iluminada devido à passagem 
de um obturador solidário com o eixo do modulador. 

Um ensaio da estabilidade do movimento do mo- 
dulador, feito com o auxílio deste sistema de medição 
do período, permitiu concluir que as variações de velo- 
cidade de um ciclo para outro eram desprezáveis. 


2.3 — Detecção e registo da modulação neutrónica 


A detecção da modulação neutrónica é feita por 
meio de uma câmara de ionização não compensada do 
tipo Westinghouse WL 6937. 

À corrente fornecida pela câmara dá entrada num 
amplificador corrente-tensão funcionando como adapta- 
dor (fig. 2), à saída do qual se procede à compensação 
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AS artar Se elevada impedancia 


AZ,ÃI] ampiitradores 


Fig. 2 -- Esquema do canal de detecção de neutrões 


quase total da componente contínua do sinal por meio 
de uma contra-tensão apropriada, O sinal é então ampli- 
ficado e digitalizado por intermédio de um conversor 
tensão-amplitude (o conversor de um selector de ampli- 
tudes de 400 canais, Intertechnique SA 40 B). 

O registo da modulação é feito através de uma 
montagem electrónica, esquematizada na fig. 3, que 
permite: 
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Fig. 3 — Esquema do equipamento de registo da modulação 
neutrónica 


— efectuar uma amostragem instantânea do sinal 
com um período de amostragem escolhido: 
(1-2-4-8 ) =<(10-5— 10—4— 103 10-2 101 

— 100 — 101 — 102) 5; 


— armazenar as sucessivas amostras nos 400 ca- 
nais da memória do selector (os impulsos do 
relógio que determinam os instantes da amos- 
tragem comandam também o avanço do endereço 
dos canais); 


— repetir a análise o número desejado de vezes, 
somando as informações provenientes de cada 
uma, nos canais correspondentes (para tanto, o 
início de cada período de análise é comandado 
por um impulso de sincronização fornecido por 
um dispositivo de referenciação solidário com o 
eixo do modulador). 


Este método de registo permite melhorar a preci- 
são dos resultados, o que é particularmente relevante 
atendendo a que a amplitude da modulação é fraca. 


3 — REALIZAÇÃO EXPERIMENTAL 


3.1 — Configuração 


O estudo da função de transferência do RPI incidiu 
sobre a configuração representada na fig, 1. 

O núcleo do reactor é constituído por elementos 
de combustível do tipo MTR, com 12 placas curvas de 
espessura igual a 1,78 mm (0,070”), formadas por uma 
liga de urânio e alumínio (45% de urânio, em peso) e 
com bainhas de alumínio. O urâneo é enriquecido a 20% 
no isótopo 235. 

O reactor pode funcionar com o núcleo imerso 
em qualquer das duas secções da piscina mas é mais 
frequente utilizar a menor porque nela estão situados 
todos os dispositivos de irradiação (designadamente, os 
tubos de irradiação e a coluna térmica). Nesta secção, 
onde foi feito o estudo, uma das faces do núcleo con- 
tacta com a coluna térmica. 
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O modulador ocupava uma posição quase central 
e o detector de neutrões estava instalado sobre a gre- 
lha (fig, 1). 


3.2 — Medição 


A determinação da função de transferência por 
modulação da reactividade implica: 


a) A estabilização do reactor com o modulador a 
rodar a uma dada velocidade. 


b) O registo da modulação neutrónica logo que o 
regime de equilíbrio se estabelece. A memória 
do selector de amplitudes contém, ao fim de 
R períodos de análise, uma imagem discreta da 
resposta do reactor à excitação introduzida (o 
canal i registará uma contagem n; queéa 
soma das amostras correspondentes obtidas em 
cada um dos períodos). Esta imagem é perfu- 
rada numa fita de papel por intermédio de um 
«chassis» de comando Intertechnique BK 21 e 
de uma perfuradora Tally 420. 


c) A medição do período do movimento (8 2,2). 


d) A introdução dos valores experimentais num 
programa de cálculo que fornece, para cada 
frequência estudada, a amplitude relativa e a 
fase da função de transferência do reactor. 


As medidas foram feitas em vários dias a potên- 
cias da ordem de 50 W. 


4 — TRATAMENTO DOS DADOS 
4.1 — Análise de Fourier da resposta do reactor 


A imagem da modulação neutrónica é uma tensão 
v It) que se apresenta sob a forma de uma sucessão de 
amostras v(t, ) colhidas em instantes separados por 
intervalos de tempo iguais. Como resultado desta amos- 
tragem instantânea aplicada à função contínua v (t), 
obtém-se, em teoria, uma nova função [10]: 


+ + 4 ê 
(= E AtvidÃoO St—iAt) 


+ oo 
Atvit) E SMt-iAt) 


— 0 


viti d (t) 


(4.1) 


com 


d. (1) = At 's 


St AU. (4.2) 


As amostras v (i At) são taxas de contagem «instantã- 
neas», isto é taxas de contagem medidas em intervalos 
de tempo muito curtos comparados comAt. Isto signi- 
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fica que se dispõe de valores de v (i A t) multiplicados 
por um factor de escala e que a introdução do factor A t 
na expressão (4.1) só se justifica por uma questão de 
comodidade. 

A transformada de Fourier de v+ (t) vem dada por 


+ a 
Frio +) — quad 
vt(g=Efv (D)=V(D'D (= E Vif*ME Pri 


— qo 


+ mo += 
=IVt—-S5J=vV(n+ 3|VE-SDS+AVEA O |. 
dá (1) 4 | v AO + + 


(4.3) 


Se o fenómeno estudado não apresentar (ou for despre- 
zável para) frequências superiores a uma frequência má- 
xima f,, e se o período de amostragem for escolhido 
de maneira que 


f.>2f4 custa Pr=(At) Stu, (4.4) 
tem-se; 
VH(f =v(f) 


No caso presente, o sinal é periódico e a sua de- 
composição em série de Fourier fornece o seguinte 
resultado: 


v(b=2|A,|eJ%h elZ=ht/T 
h | | (4.5) 
cuja transformada de Fourier é: 
V(fj= 2 [A | ed 5a ai 
dE E (4.6) 
h h r 


Por conseguinte, tendo em conta (4.4), pode escrever-se: 


At (=A, (f) (4.7) 

Por outro lado, se as frequências dos harmónicos 
de v(t) são inferiores a (2 At) —1, não se observa qual- 
quer efeito de dobragem e obtém-se a fase da função de 
transferência também sem qualquer correcção. 

Basicamente, o princípio do tratamento dos dados 
relativos à modulação neutrónica consiste em efectuar 
uma decomposição em série de Fourier da porção de 
cada registo correspondente ao período de excitação 
previamente determinado (8 2.2). 

O conteúdo de cada canal da memória do selec- 
tor, ao fim de R períodos de análise (fig. 4), é 


Y : a R 
n(D) =, ER DO ron) 


(4.8) 


Calcula-se o maior inteiro contido em T/ At (seja N) e 
efctua-se a decomposição de Fourier sobre as N 
amostras: 
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ntt;) 


PE np 2 a 


- —— a —— + 


Fig. 4 —— Amostragem do sinal 


2 ; 
h h RT ias 1 ri 1 
ama Zn; ej2rhiAt'T (4.9) 
| 
A amplitude do harmónico é dada por 
12 : 
| ” 4.10 
e Oo seu argumento por 
| 
Mp =arc tg S,/C, | (4.11) 


Note-se que vp é a fase da resposta do reactor em rela- 
ção à referência imposta pelo impulso de disparo dos cl- 
clos de análise (8 2.3). 

O nível médio do sinal vem dado por: 


Dri (4.12) 


NEÇEr 


- ta Z 


e permite calcular a amplitude relativa do harmónico h: 


a, = Ap'vo (4.13) 


Um programa escrito em linguagem Fortran para 
um computador PDP 15/20 permite calcular th Cho Sp 
A, ?n € à, para os três primeiros harmónicos, a par- 
tir dos valores de T, R, V;, K en,. 
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4.2 — Determinação da função de transferência 


Para determinar a função de transferência de um 
“reactor não é imprescindível excitá-lo de uma maneira 
rigorosamente sinusoidal, o que, de resto, é pratica- 
mente impossível. Basta introduzir uma excitação perió- 
dica (de fraca amplitude) que pode ser considerada 
como uma soma de sinusoides puras, cada uma das 
quais contribui para a resposta do reactor. Às compo- 
nentes fundamentais obtidas por análise de Fourier da 
resposta e da excitação permitem calcular a função de 
transferência. 

Acontece porém que não é muito importante deter- 
minar o ganho da função de transferência em valor 
absoluto, o que permite passar sem o conhecimento 
da forma exacta da excitação (apenas há que garantir 
que a amplitude é fraca). Nestas condições, basta consi- 
derar, para cada frequência, o fundamental da decom- 
posição da resposta; concretamente os valores de a, 
e v, para h= 1. 

Para obter os valores absolutos da fase da função 
de transferência do reactor, importa ainda determinar 
a diferença de fase entre a referência que marca o iní- 
cio dos ciclos de análise e a excitação. Para isso, faz-se 
rodar o modulador nos dois sentidos, à mema frequên- 
cia: a semi-diferença dos valores de vw assim medidos 
permite calcular esta diferença de fase; por seu turno, 
a semi-soma daqueles ângulos fornece directamente a 
fase da função de transferência correspondente à fre- 
quência estudada. 

A partir da curva experimental da amplitude rela- 
tiva da função de transferência, medida com o reactor 
crítico, é possível determinar a frequência de corte, f.. 
e a partir desta o valor de x 


Também é possível representar os valores de tg 
em função de o —> »y (bm 6 0 maior 2 ): obtém- 
-Se uma recta cuja inclinação é «(1 


4.3 — Correcções 


4.3.1 — Influência da banda passante do ampli- 
ficador 


O amplificador associado ao detector de neutrões 
pode deformar os sinais em amplitude e em fase, Sendo 
F. a frequência de corte do amplificador, há pois que 
efectuar, em princípio, as seguintes correcções: 


Sobre a fase: = tactg(f/F,) 


Sobre a amplitude: A, Ma | t+ (F/ Fo)? | H2 
(v e m significam respectivamente verdadeiro e medido). 
Note-se que, para f=F./10, a diferença de fase 
introduzida pelo amplificador é ainda de — 6º enquanto 
que a razão entre as amplitudes é praticamente igual a 1. 
No presente trabalho tem-se F. — 1,6 kHz pelo 
que aquelas correcções são praticamente desprezáveis. 
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4.3.2 — Influência do tipo de registo 


Como se viu no 8 4.1, a amostragem instantânea 
à qual é submetido o sinal vít) não provoca qualquer 
deformação na fase e amplitude respectivas, desde que 
o espectro do sinal possa ser considerado limitado e 
desde que se escolha uma frequência de amostragem 
tal que fwy > fm: 

No caso presente, tem-se sempre 


200 At < T< 400 At, 


200 100 
ou seja E > fy > T 

Por conseguinte, a frequência de amostragem é 
sempre nitidamente superior à frequência dos harmóni- 
cos do sinal, os quais, de resto, são praticamente des- 
prezáveis a partir da 3.º ordem (numa decomposição 
típica correspondente a uma frequência do fundamen- 
tal igual a 0,035 Hz, a amplitude do harmónico 2 vale 
3,2% da amplitude do fundamental, ao passo que a 
amplitude do harmónico 3 só representa 0,9%). 

Nestas circunstâncias, os valores da amplitude e 
da fase do fundamental de v(t) não foram sujeitos a 
qualquer correcção proveniente do tipo de registo. 


5 — PRECISÃO 
5.1 — Generalidades 


As contagens n, revelam flutuações devidas à 
natureza estatística 


— das reacções nucleares que ocorrem no reactor, 

— dos fenómenos ligados à detecção dos neutrões, 

—do ruído electrónico do canal de detecção, 

— da conversão analógico — digital, 

— da reprodutibilidade das condições de amos- 
tragem. 


Admitindo que as quantidades n ,, são variáveis gaus- 
sianas, as distribuições de probabilidade de C, e sS h 
são também gaussianas, dado que estes termos do desen- 
volvimento de vit) são obtidos através de uma série de 
operações lineares que incidem sobre os valores de n,, 
Nesta hipótese, os desvios-padrão de C, e S,-9cecç- 


podem ser considerados como as incertezas que afec- 


tam o valores respectivos. 

Além de estimar a precisão com que os resultados 
de uma experiência representam as esperanças matemá- 
ticas de C, e S, , importa determinar os erros devidos 
ao tipo de análise a que é submetido o sinal. 


5.2 — Erros devidos às flutuações neutrónicas 


Pretende-se determinar a precisão das medidas 
da amplitude e do argumento do fundamental de v(t), 
examinando-se a influência das flutuações estatísticas 
do sinal nas seguintes condições: 
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— ausência de modulação; 

— análise contínua (a decomposição de Fourier 
incide sobre um sinal ainda não submetido a 
qualquer amostragem). 


De acordo com exposto na ref. [9], conclui-se que: 


a) As medidas C e S representam as respectivas 
esperanças matemáticas a mais ou menos 
2ºc = 2º, com uma probabilidade de 95%. 

b) O desvio-padrão da medida da amplitude do 
fundamenta! vem dado por: 


| P(1 T) | 12 
G — Gm» = [ n = 
Â e 5 RT 


c) O desvio-padrão da medida do argumento do 
fundamental é 


(5.1) 


Ejs mico (5.2) 


d) A razão sinal/ruído, na hipótese de se despre- 
zar o ruído da aparelhagem (o que no caso 
presente se revelou perfeitamente razoável ex- 
cepto para f = 5O Hz), pode exprimir-se sob 
a forma 


A | 
— = Ac |H(h| , dida 21 * 
A Cm tsHo|S 


| 


(5.3.) 


Facilmente se conclui que para melhorar esta razão há 
interesse em utilizar um detector com a maior eficiên- 
cia possível e em aumentar o tempo de registo e/ou 
a potência do reactor (sobretudo quando a eficiência é 
pequena e se trabalha a frequências elevadas pois, então 
para registos de duração igual, S/R diminui de 6 
dB/oitava na banda de frequências em que H(f) varia 
desta maneira). 


5.3 — Erros devidos ao tipo de análise 


No 8 5.2, examinou-se a influência das flutuações 
estatísticas do sinal sobre a precisão dos valores me- 
didos, supondo que a análise era de tipo contínuo. Na 
realidade, o sinal v(t) que sai do amplificador corrente- 
-tensão é submetido a uma amostragem periódica, sendo 
cada amostra o resultado de uma conversão analógico- 
-digital. Finalmente a decomposição de Fourier vai inci- 
dir sobre a representação discreta de vít) assim obtida. 
Importa, por conseguinte, analisar as implicações deste 
tipo de análise na precisão dos resultados. 


a) Conversão analógico-digital 


A conversão analógico-digital acarreta uma quan- 
tificação do sinal: os valores contínuos de vit) vêm 
representados por valores discretos nít,). O conversor 
impõe um limite superior a n(t,) que depende, em 
particular, do coeficiente de conversão escolhido — e 
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cujo inverso é afinal a unidade da escala de medição. 
O erro relativo introduzido pela quantificação é da 
ordem desta unidade e no caso presente, é, em geral, 
inferior a 1%, 


b) Condições de amostragem 

As deficiências na reproductibilidade das condições 
de amostragem podem ser imputadas: 

— ao sistema de disparo dos ciclos de análise; 

— ao relógio que comanda os instantes de amos- 

tragem. 

As características de um e de outro permitem toda- 
via desprezar erros provenientes de tais origens. Em 
particular, verificou-se com um relógio de quartzo que 
a estabilidade dos valores dos períodos de amostra- 
gem era excelente. 

Por outro lado, a soma dos nít,) correspondentes 
aos diversos ciclos de análise é perfeita e não impõe 
restrições. 


c) Decomposição de Fourier 


A decomposição de Fourier do sinal incide sobre 
mais de 200 amostras por período de modulação, admi- 
tindo-se que daí não resulta qualquer influência sobre 
a precisão dos resultados (5 4.3.2). 


5.4 — Cálculo de erros 


A expressão (5.3) permite calcular o erro relativo 
do valor A = A, , expressão (4.10), da amplitude co 
fundamental para condições de medição determinadas, 
isto é, para valores conhecidos de F,:, 4: Re f= 1/T. 
Recorda-se que se tem y = 0,8, 

E=32: 10" p (5.4) 
e que o número de acontecimentos úteis registados pela 
câmara pode ser estimado [7] por meio de 


1,6. 10-19 E, 1 
No = i DE F 
e Wn) = 


(5.5) 


em que E, — 1,47 MeV, E, — 26 eV e F, — 0,6. 


O desvio padrão de A vem então dado por 


R 
ma MAM. Sd) (5.6) 


. 


Para obter o desvio padrão da amplitude relativa 
a — A/v, há que ter em conta, finalmente, a impreci- 
são de v, . 

O valor de R/S fornece directamente Sw e para 
determinar o erro da fase ? da função de transferência 
basta introduzir a incerteza relativa à diferença de fase 
entre a referência que marca o início dos ciclos de aná- 
lise e a excitação (8 4.2). 
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6 — RESULTADOS 


Os valores experimentais da amplitude relativa e 
da fase do fundamental da resposta do RPI, obtidos a 
baixa potência na configuração estudada, estão repre- 
sentados nas figs. 5 e 6. Os erros relativos indicados 
na fig. 5 e as barras de erro representadas na fig. 6 
corresponcem a | 2965,/a ea + 26, respectiva- 


e CRVA ERPERIME TAL 
[valores normalizados A=" para f=1m8) 
| ——— MODELO PONTUAL , 
| (B, = 75010 A 2810 3) 


Í 
o. + 4 Am e Deo Que Apr am mm q 44 Ls s4 


mente, e são típicos da gama de frequências em que se 
situam. A sua determinação foi feita nas condições ex- 
postas no capítulo 5 e constata-se que 25,/a < 0,5% e 
2c0cp < 1,0º em quase toda a banda de frequência estu- 
dada. Verifica-se ainda que, para frequências superiores à 
frequência de corte do reactor, Os erros aumentam à 
medida que H(f) diminui. 


E rreT ' pq I E RES E em) na sas pude é (EE SS potó CR (O (CU I aa 
| 
4 


10 5 (0) f ima) + 10 


Fig. 5 — Amplitude relativa da função de transferência do RPI à potência zero 


60 — 
—— Ciaya EXPERIMENTAL 
mise MODELO PONTUAL 
E] ç 
(Re =780 0 ; A =SEN 5) 


Bo : 


-— e--—s 


Fig. 6 — Fase da função de transferência do RPI à potência zero 
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QUADRO | 


Grupo T, (s) A; (s—1) B, /B 
1 54,51 + 0,94 0,0127 + 0,0002 0,038 + 0,003 
2 21,84 + 0,54 0,0317 + 0,0008 0,213 + 0,005 
3 6,00 + 0,17 0,115 + 0,003 0,188 + 0,016 
4 2,23 + 0,06 0,311 + 0,008 0,407 -+ 0,007 
5 0,496 + 0,029 1,40 + 0,081 0,128 + 0,008 
6 0,179 + 0,017 3,87 + 0,369 0,026 -+ 0,003 


Aos valores experimentais foi ajustada uma curva 
desenhada a partir dos valores calculados em teoria 
pontual utilizando os parâmetros característicos da emis- 
são dos neutrões retardados que se encontram coligidos 
no Quadro | [81]. 

Da comparação dos valores experimentais com os 
calculados infere-se que 


— O tempo de geração dos neutrões instantâneos 
corespondente à curva que melhor se ajusta aos 
resultados experimentais é 


1=(56+2). 107º &: 


a = 135,7 s SEÉ + 0,5%) e o valor estimado 
para $, é (760 + 23). 10 *; 


— para frequências superiores à frequência de corte 
do reactor, 


f. = 21,6 Hz (+ 0,5%), 


o modelo pontual não é suficiente para repre- 

sentar o comportamento dinâmico de um reactor 

nuclear pois os efeitos de espaço deixam de ser 
desprezáveis [9]. 

Uma experiência efectuada segundo o método 

Rossi- x, numa configuração mais pequena e sem o 

elemento central retirado (actual posição do modulador), 


fornecera o resultado [5]: A= (484-4).10 ss. 

A determinação de A pelo método exposto no 
5 4.2 revela-se delicada devido à presença dos harmó- 
nicos espaciais. 


CONCLUSÕES 


O método de modulação da reactividade permite 
determinar o ganho e a fase da função de transferên- 
cia de um reactor, no ponto em que se efectua a me- 
dição. A partir daqueles valores, é possível deduzir o 
tempo de geração dos neutrões instantâneos desde 
que se considere conhecida a fracção efectiva dos neu- 
trões retardados 5. 
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Com o modulador concebido para esta experiên- 
cia foi possível estudar, com uma boa estabilidade, uma 


banda de frequências compreendida entre 5.10 —3 Hz e 
65 Hz. A técnica experimental adoptada para a detec- 
ção e análise da resposta do reactor revelou-se muito 
eficaz. E os resultados obtidos para a função de transfe- 
rência do RPI, à potência zero, forneceram um valor 
para à = (56+2). 109 s com 5, = (760+423).10 * 
Constatou-se mais uma vez [9] a influência dos har- 
mónicos espaciais no comportamento dinâmico dos 
reactores, a qual não é tida em conta pelo modelo 
pontual usualmente utilizado nos cálculos teóricos. 
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SIGNIFICADO DOS SIMBOLOS 


A, — amplitudo do harmónico h do desenvolvi- 
mento de n,, [t) em série de Fourier. 

a, — amplitude relativa do harmónico h. 

C, (ti— concentração dos precursores de neutrões 
retardados do grupo í no instante t, 

C, — termo em coseno correspondente ao harmó- 
nico h do desenvolvimento de n, em séria 
de Fourier. 

E, — energia necessária para produzir um par de 
jões. 

E, — anergia de uma partícula z emitida numa 
câmara de ionização com depósito de boro. 

F — número de cisões ocorridas no reactor por 
unidade de tempo. 

F, — fracção média da energia da partícula , que 
provoca a ionização do gás da câmara. 

F. — frequência de corte do amplificador. 

F(fit)] — transformada de Fourier da função fit). 
É -1/T frequência, 


Ea (At)! frequência de amostragem. 


f. — frequência de corte do reactor. 


f, — frequência do harmónico h. 
ty (2At)”! frequência de Nyquist. 
f,,— frequência máxima do espectro de vit). 
| Hj) | — ganho da função da transferência do reactor. 
Hljrm) — função de transferência do reactor. 
H, — H(jhw) = H,. exp (jt,). 
| — corrente fornecida pela câmara de ionização. 
jP=—1, 
K — coeficiente de conversão do analisador. 
L(f(t) ; — transformada de Laplace da função f(t), 
| — vida dos neutrões instantâneos. 
m — número de grupos de neutrões retardados, 
Nilo) -Lt(nit)) 
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N T/At. 
n(t) — densidade neutrónica no instante t, 
ny, -mljho) =n,. exp (jz,). 
n —n,lt,) conteúdo do canal i ao fim de R pe- 
riodos de análise. 
no Censidade neutrônica de equilibrio. 
P(f) — densidade espectral energética de um ruído, 
R — número de periodos de análise. 
Rijo) = fe (t)) 
S/R — razão sinal/ruido. 
Sh — termo em seno correspondente ao harmónico 
h do desenvolvimento de n,. (t) em série de 
Fourier. 
T — período fundamental da modulação. 
Tá — período de análise. 
IF, — semi-vida dos precursores de neutrões retar- 
dados do grupo f. 
T; — 400 À t duração do registo. 
t — tempo, 
É = instante correspondente ao canal i. 
Vit) — F(fvít)) 
V., — contra-tensão. 
vit) — tensão-imagem da modulação neutrónica. 
vt (t) — função obtida por amostragem de vit). 


a = BIA 

C a -: 

De — valor cficaz da fracção total dos neutrões 
retardados. 


3. — fracção dos neutrões retardados do grupo 1. 
ER aficiância dos neutrões retardados do grupo i, 
An it; — variação de n(t) de fraca amplitude. 
Ae lt; — variação de clt) de fraca amplitude, 
At — período de amostragem. 
5 [ti — função delta de Dirac. 
: — eficiência do detector. 
+ — tempo de geração dos neutrões instantâneos. 
“;— constante de desintegração dos precursores 
de neutrões retardados do grupo i. 
“ — número de neutrões emitidos por cisão, 
: |t| — reactividade no instante t. 
en po exp (iv). 
*, — reactividade no instante t= 0. 
"« — desvio padrão de x, 
b — fase da função de transferência do reactor. 
", — argumento de H, 
“y — argumento de nh 
p=<up—1)>peu>? 
“n — argumento de « 
w— 2 =f pulsação. 
- — média de um conjunto. 
* — produto de convolução. 


h 
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SOME SPHERICAL-SIMEIRY BOUNDARY-VALUE PROBLEMS IN STEADY- 
-STATE ANISOTHERMAL ELASTIC MEDIUM 


RESUMO 


O objectivo deste artigo é estudar uma classe de 
problemas inter-relacionados, respeitantes à determina- 
ção de temperatura, vector de fluxo de calor, vector de 
deslocamento, tensor de deformação, tensor de 
tensão e vector de tensão, de um meio linear, homogé- 
neo, isótropo e elástico, em regime permanente, com 
condução de calor, para 4 geometrias esféricas — esfera 
oca, camada esférica, esfera maciça e cavidade 
esférica num meio infinito — , e condições na fronteira 
especificando a deslocação superficial, a pressão ou a 
tensão tanjencial, e a temperatura, o fluxo de calor ou 
condições de convecção. 

Alguns problemas de simetria esférica às condi- 
ções na fronteira em regime permanente num meio 
anisotérmico e elástico. 


NOTATION 
E = scalar density Young's Modulus 
q — scalar Poisson Ratio. 
a = scalar density Thermal Expansivity. 
C,. C, = scalar densities Specific Heat at 


Constant Pressure and Volume, 
k = scalar density Thermal Conducti- 
vity. 
à jo m=1,2,3= running indices; o, = dt 
); = Kronecker Delta multiplicity (com- 
ponents are 1ifi=j, and Oifi ij). 
D, = Covariant Derivative. 


= contravariant vector Displacement 
Vector. 


S; = mixed tensor of order 2 Strain 
Tensor. 


T; = mixed tensor density of order two 
Stress Tensor. 

T,= covariant vector density Stress 
Vector. 

T = scalar Temperature (scale zero at 
zero-Stress-Strain). 

Ji = contravariant vector density Heat 
Flux. 

n, = covariant vector Normal to a regu- 

lar surface, 
r, 9, v = spherical coordinates. 
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LUIS MANUEL BRAGA DA COSTA CAMPOS 


SINOPSYS 


The purpose of this article is to study a class of 
interrelated problems, concerning the determination of 
the temperature, heat flux vector, displacement vector, 
strain tensor, stress tensor and stress vector, of a linear, 
homogeneous, isotropic elastic medium, in steady-state, 
with heat conduction, for 4 spherical geometries — 
hollow sphere, spherical layer, solid sphere and spheri- 
cal cavity in infinite medium — , and a range of boun- 
dary conditions, — given surface displacement, pressure 
or tangential stress, and temperature, heat flux, or con- 
vection conditions. 


u = radial physical component of dis- 
placement vector. 
j = radial physical component of theat 
flux. 
Tis Tha To = spherical physical components of 
Stress Vector. 


Too Te, Tec, Tom Ty Toy = spherical physical compo- 
nents of Stress Tensor, 
S rr 370, Sr, So Spy Syy = spherical physical compo- 


nents of Strain Tensor. 
A, B, C, D = arbitrary constants of integration. 
R, R,, R; — boundary surface's radiuses. 
To» T,.T ; = boundary surface's temperatures. 
Jor Jj,s Jo — boundary surface radial physical 
component of Heat-Flux, 
Pp, Pj, Pp = boundary surface pressure. 
t, tj, !, = boundary surface tangential stress. 
Uo, Uj, Us = boundary surface radial physical 
component of displacement. 


| — THE FOUNDAMENTAL EQUATIONS 

1.1 — The foundamental equations in tensor from 
Considering a linear, homogeneous, isotropic, ani- 

sothermal elastic medium, in steady state, in the 


absence of gravitic or eletromagnetic fields, and chemi- 
cal reaction or phase change, of scalar density Young's 
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Modulns E, scalar Poisson Ratio q, scalar density Ther- 
mal Expansivity a, scalar densities Specific Heats at 
Constant Pressure and Volume, C, and C,, and scalar 
density Thermal Conductivity k, the following basic 
equations hold: 

(Il) Relation of the contravariant vector density 
Heat Flux, j ', to the scalar Temperature, T, — Fou- 
rior's Law 


J' -k D*T (1) 


(ll) Defenition of the mixed tensor of order two, 
Strain Tensor, S; from the contravariant vector Displa- 
cement Vector, u ; 


S;= (1/2) (D;u'+D'u,) (2) 


(Ill) Relation between the covariant vector density 
Stress Vector, T, , and the mixed tensor density of 


order two Stress Tensor, T; at an infinitesima| regular 
element of area of covariant vector Normal, n,, 


T.=-Tin. (3) 


i a 


(Iv) Mixed tensor density of order two Stress 
Tensor, LÉ , as function of the mixed tensor of order 


two Strain Tensor, Ss] and the scalar Temperature, 
— Hooke's Law, 


j=(E/MM1 —q)) S;+ (Eq/((1 + q) (1-2q)) SMS; + 
+ (aE/(3 (1-2q))T (4) 


(V) Second-order vector and scalar linear homo- 
geneous differential equations, with constant coeficients, 
satisfied by the contravariant vector Displacement Vec- 
tor, u!, and scalar temperature, T — momentum and 
energy equations 


((1-29)/(1+q)) D Diu! + (1/(1+q)) D'D u! = 
=(2/3)a D'T (5) 


C,dAT+UC,— Cla) dD;u! =k D,D'T (6) 


|.2 — Foundamental Equations with Spherical-Simetry 


Assuming spherical-simetry, all components of the 
displacement vector vanish, except for the radial com- 
ponent that depends only on r; the temperature also 
depends only on r, 

Rewriting the preceeding equations, (1) to (6), in 
physical components, in spherical coordinates, and using 
the above stated assumptions, we conclude that: 

In a (1) linear, homogeneous, isotropic, anisother- 
mal elastic medium in steady state, in the absence of 
gravitic or eletromagnetic fields, chemical reaction or 
phase change, of (2) Young's Modulus E, Poisson 
Ratio q, Thermal Expansivity a, and Thermal Conducti- 
vity k, the following relations hold 
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(|) Components of the Heat Flux vector 
JS ()=J(rn)=—k dT/dr (7) 


Jo (=), ()=0 (8) 


(ll) Components of Strain Tensor as function of 
Displarcement 


S, (r) = du/dr (9) 
Sqo lr) = Sou lr) = u/r (10) 
Sot)=S yl)=S (r)=0 (11) 


(Ill) Relation between Stress Tensor and Vector, 
for spherical surface element 


T, tr) = Tr) (12) 
To lr) = Too dr) (13) 
To lr) =Too tr) (14) 


(IV) Generalized Hooke's Law 


Te = (EM = q)/1 + ql —29))) S, + (Eq/((1 + 
+ q)(1 — 2q))) (Sg + Sos) + (E /(3(1 — 2q))T 


To, = (EM — q)/U1 + att — 29))) Sgg + (Ea/((9 
qM1 — 29))) (8, + 84) + (aE/(301 — 29))T 


Teo = EM — q)/1 + qt — 2q))) Soo + (Eg/((1 
qM(1 — 29))) (8, — Syg) + (aE/(3(1 — 2q))T 
(15) - (16) - (17) 


Tio = (EMI +) S, (19) 
Teo = (EMI +) S,, (18) 
Too = (EMI + q)) Sy (20) 


(V) Momentum and Energy Equations 


(3(1 —- q)/(1+q)) d/dr ( (1/r2)d/drr2 ulr) ) = a dT/dr 
(21) 


r2 doT/dr2 + 2 rdT/dr4t2 T=0 (22) 


|| — GENERAL SOLUTION 
1 — General Integral of (21) - (22) 


We shall now obtain the general integral of the 
system of two linear homogeneous diferential equations 
of order 2, dependent variables u and T, independent 
variable r, formed by the equations of momentum and 
energy, (21) - (22). 
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We remark that equation (22), the Equation of 
Energy, is a linear homogeneous diferential equation of 
order two, with a single dependent variable T, and is 
of form known as Euler-type, witch the change of inde- 
pendent variable 


s=logr r=es (23) 


reduces to a ljnear homogengous diferential equation of 
order two, with constant coeficients 


d2T/ds2 + dT/ds—-2 T=O0 (24) 


whose General Integral can be expressed as a linear 
combination of two linearly independent Particular 
Integrals: 


T(IS)=A T,(s)+B T.(s) (25) 


with A and B two arbitrary constants of integration. 
Remarking that equation (24) can be wiriten 


( d/ds-1) ( d/ds+ 2 ) Tiís) =0 (26) 


two particular solutions are 
(d/ds-1) T(s) = O 


(d/ds+ 2) Tí(s)=0 


(27) - (28) 


to with correspond the two linearly independent Parti- 
cular Integrals 


T, (s)b=es Ta (s)=0e—2 (29) 


Replacing (29) in 25), and using (23), we con- 
clude that the general integral of the equation of energy, 
(22), i. e., the Temperature, is 


T(r) = A r+B r-—2 (30) 


in witch A and B are arbitrary constants of integration. 
Replacing (30) in (21), we obtain the second- 
-order linear inhomogeneous diferential equation satis- 
fied by the displacement 
d/dr r-2 d/dr r2 ulr) = (a(1 + q)/(3(1 — q)) 
(A —- 2 Br-3) (31) 


or, integrationg once, and introducing the arbitrary 


constant of integration C 


d/dr r? ulr) = (a(1 + q)/(3(1 — q)) 
(Ar? +B) 4-3C r2 (32) 
or, integrating once more, and introdueing the arbi- 


trary constant of integration D, the expression of the 
Displacement is 


u(r) = (a(1 + q)/(3(1 — q)) 
(Ar2/4+Br-1)) +Cr+D q-—2 (33) 
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[|.2 — General Solution 


Hence, a (1) linear, homogeneous, isotropic, ani- 
sothermal elastic medium, in steady state, in the 
absence of gravitic or eletromagnetic fields, chemical 
reaction or phase transition, (2) of Young's Modulus E, 
Poisson Ratjo q, Thermal Expansivity a, Thermal Con- 
ductivity k, (3) with physical spherical simetry, satis- 
fies the following equations, in witch A, B, C and D 
are arbitrary constants of integration, specified by the 
boundary conditions. 


(|) Temperature and Displacement, from (30) and 
(33) 


T(r)=A r+B r-2 (34) 


u(r) = (a(1+q)/(12(1—-q)) Ar2+Cr+(a(1+ 


+q)/(9(1 —q)) Br-l+D r—2 (35) 

(ll) Heat Flux spherjcal-physical componentes, 
from (7) - (8) and (30) 

J(r)=) ()=—-k A+Zk B r—3 (36) 

Jo ()=J; ()=0 (37) 


(Ill) Strain Tensor spherical-physical components, 
from (9) - (106) - (11) and (35) 


Sall=(alt+a/(6(1-q) Ar+C-(a(1+ 
“a(s —q))Br-2-2D r-3 (38) 


Soo (r) = Say lr) = (a(1 + q)/(2(1 — q)) Ar+ Ch+ 
(a(1 + q)/(3(1 — q)) Br—-2+ Dr—3 


3.6 (1) = 3,9 (r) = S by (r) = 0 (39) - (40) 
(IV) Stress Tensor and Vector (for spherical sur- 


face elements) spherical-physical companentes, from 
(12) to (20) and (38) to (40); 


El — q) E 
T —— Ba) É .——e—e—e—e—e—e 0 PAO fo mem 
A od dd "fa 
2a Eq 2E 
ES mare irreeça D m A  cem, Ai RD, UM) 
> [4 =0) Ti — 4) 1+q 


Told=To (WM=To(M=Tro (= 


aE (5 —2q) E 
=" A+ -—— € 
12 (17 — 9) (1 — 29) 1+aq 
2a E E 
+ - —— Br-: + — Dr? 
do gq ii 
Tool) =Try ()=To ()=0 (43) 
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1.3 — Boundary conditions 


The constants of integration A, B, €C and D, in 
Formulae (34) to (43) are determined by the boundary 
conditions, i. e., specified by the conditions satisfied 
by the temperature and displacement, at given medium 
limiting or boundary surfaces. 

Whatever the geometric shape of the medium, in 
order to preserve spherical - simetry, the boundary sur- 
faces must be concentric spherical surfaces. We shall 
consider 4 geometries: 

(1) Hollow Sphere of inner and outer radiuses R, 
and Rs. 

(Il) Spherical Layer of radius R and thickress e, 
e < < R; equivalent to (i) with 


R=(R, +R,)/2 e=R,-R, (44) 


= 
(ll) Solid Sphere of radius R. Equivalent to (1) 
with R|, = O and R, =R. 
(Iv) Spherical Cavity of radius R in infinite 
medium equivalent to (1) with 


R=R, R; — o (45) 


To preserve spherical-simetry, the physical varia- 
bles must be constant at these boundary surfaces. As 
example of the most common boundary canditions, at 
one sutch surface, of radius R , , we give the following, 
we shall use later: for Temperature 

(|) Specified Temperature T, at boundary surface 

TQARD)=T, (46) 


(H) Heat Flux J, through the unit are of boundary 
surface in unit time, in the direction of the normal 


diád=a, (47) 


(Il) Convection of heat through the boundary sur- 
face, from a fluid at free temperature T.,, 


CIMR,)=h4TIR|)-TA) (48) 
h being the Surface Thermal Conductivity. 
for the Displacement 
(Il) Displacement u at boundary surface 
u(Rj)=u, (49) 
(ll) Pressure p at boundary surface 
T, (R,)) =p, (50) 
(Il) Tangential Stress T, at boundary-surface 
TA(RG)I=TIR)=?, (51) 
Combining these boundary conditions, we would 


obtain 3 x 3 = 9 distinct sets of boundary conditions for 
each boundary surface, thus 9 X 9 = 81 distinct pro- 
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blems for geometries (Il) and (ll), hollow sphere and 
spherical layer, and 9 distinct problems for geometries 
(Il) and (IV), solid sphere and spherical cavity, although 
we shall only give 3 + 1 examples of the former and 
3 + 3 examples of thelatter. 


ll — SOLID SPHERE 
HH.1 — Problem 1 


We consider a solid sphere of radius R, at pressure 
p and temperature T, at its surface. 

For all solid sphers problems, since the tempe- 
rature and displacement must be everywhere finite in 
the sphere, namely at its center, r= O, from (34) and 
(35) it follows that 


B=0=D (52) 

The boundary conditions for this problem are 
Te AR)=p TIR)=To (53) - (54) 
witch, replacing (41) and (34), using (52), and sol- 
ving for A and C, AT, /R C=(1-2q) p/E- 
— (3 q) aT, /(6(1 — q)) (55) - (56) 


From (52), (55) and (34), (36) the temperature 
and non-vanishing component of heat flux are 


TH) =TorR HMr)=—k T,/R (57) - (58) 
the former increasing linearly from O at the center of 
the spere to (54) at its surface, the latter being cons- 
tant. 


From (52) - (55) - (56) and (35), the displace- 
ment is 


utr) = ( (1-2a)p/E — ( (3 — q)/(6(1 — q)JaT,) 
r+H(1+a)/(3(1 — q)) (aT,/R) r2 (59) 


varying fromO at the center of the sphere, to 
u(R) = ((1-2q)/E) pR — aT, R/6 (60) 
at its surface, through the maximum 


Rn U3— q) q)DR-(30-q) (1— 29) /(201 + 


“ q))) (pR/aET,) (61) 
UR) =(U3—q) (1 — 29)/((8(1 + q))) pR/E — 
—(U3—-q)2/(48(1 —- q2))) aT,R (62) 


The non-vanishing components of the strain tensor, 
are, from (2) - (55) - (56) and (38) - (39 - (40) 


Se tr)= (1 +al/(6(1 — q))) aT; r/R + (1 — 
— 2q)/E)p—((3 — q)/(6(1 — q))) af, (63) 


TÉCNICA 419 


S og (r)= seo tr) = (1 + q)/(1201 — q) ate RA 1 — 
— 2q9)/E)p — (43 — q)/(6(1 — q))) af, 


The non-vanishing components of the Stress Ten- 
sor are, from (52) - (55) - (56) and (41) - (42) - (43) 


Ta t)=p+(U3-—aq)/(6(1 — q) (1 — 2q))) aET, 
(r/R — 1) (65) 


Too lr) =T )=p+ (aET,/(12(1 — q) (1 — 29))) 


eq fr 
((5 — g)r/R + 2q + 6) (66) 


The Stress-state at the center of the sphere, is 
from (65) - (66), spherical 


T, (0) =To (0) =T, (0)=p — aT, El(3 — q) /(1601 — 
— q) (1 — 2q)) (67) 


and the stress-stat at the surface of the sphere is hemi- 
sprehical specified by (53) and 


To (R)=To(R)=p—aET, /(6(1 — q) (1 — 20)) 
(68) 


1.2 — Problem 2 


We consider a solider a solid sphere of radius R, 


with displacement u, at and heat flux J, through its 
surface. The boundary conditions are 
u(R) = u, J(R) =J, (69) - (70) 


Replacing (35) in (69), and (36) in (70), using 
(52), and solving for A and C, we get 
A=-J,lk C=u,/R+((1+q) (12(1 — q))) 


as, R/k (71) - (72) 


The temperature and non-vanishing component of 
heat flux are, from (52) - (71) and (34) - (36), 


Tt) = — (J, /k) r J(tr) = J, (73) - (74) 


the former decreasing linearly from zero at the center 
of the sphere to — J, R/k at its surface. 

From (52)-(71)-(72) and (35), the displace- 
mente is 


ul) =u, /R+(1+q)/(1201 — q))) 
(aJ, /k) (R — rr (75) 


and varies from O at the center of the sphere, to (69) 
at its surface, through the maximum 


Rm = R/2+ (1201 — q)/(1 + q)) (u, k/2a9, R) 
(76) 
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Ro URm)=uU,/2R+((1+q)/(12(1 — q))) 
(a), R/2k) (77) 


From (52) -(71)-(72) and (38) - (39) - (40) the 
non-vanishing components of the Strain Tensor are 
Set) =u,/R+((1+q)/(1201-q))laJ,/k) (R— 2r) 

(78) 


Sgg (r) = Sus lr) = u,/R + (1 + q)/(12(1 — q))) (ad, 
/k) (R — f) (79) 

and from (52)-(71)-(72) and (41)- (42) - (43) the 

non-vanishing components of the Stress Tensor are 

TF, tr) =(E/(1 — 2q)Ju, /R + (aEJ, /K)((1 + q) R — 

— 2 (3 — q)r)/(12(1 — q) (1 — 29)) (80) 


Tio (r) = Too lr) = (E/(1 — 2q)Ju, /R+ aaEJ,/k) ((1+ 
+ q) R-— (5 — 2g)r)/(12(1 — q) (1 —29)) 
The Stress-state at the center of the sphere is 
spherical 
T;(0)=To (0) =T (0)= (E/(1 — 2q))u,/R + ((+ 
» q)/(12(1 — q))) aEJ, R/k (82) 


and at the surface of the sphere is hemi-spherical 


T;(R) = (E/(1 — q)lu,/R — ((5 — fa) /(12(01 — q) (1- 
— 2q))) (aEJ, R/k) (83) 


To (R)=To(R)=(E/(1 — 2q))uç,/R — ((4 — 
— Sq)/(12(1 — q) (1 — 29))) aEJ, R/k (84) 


1.3 — Problem 3 


We consider a solid sphere of radius R, tangen- 
tial stress to at its surface, and heat convection from 
a fluid at free temperature T, ; the boundary con- 
ditions are 


To (R)=To (Rj=t,  TIRI=T, (85)-(86) 


Replacing (42) in (85) and (34) - (36) in (86), 
using (82) and solving for A and C, we obtain 
A=hTo/(k + hR) C = ((1— 2q)/Et, — 

— (43 — q)/(6(1 — q))) aT, (87) - (88) 

From (87), (82) and (34) the temperature is given 
by (89), increasing linearly from O at the center of the 
sphere to (90) at its surface, and from (52), (87) and 
(36) the non-vanishing component of the heat flux is 
given by (91) 

T(r) = (hT, /(khR))r T(R) = hT, R/(k + hR) 
J(r) = — khT./(k + hR) (89) - (90) - (91) 


From (52) - (87) - (88) and (35) the displace- 
ment is 


u()=((1+aq) (1211 — q))) (ahT, /(k + hR)) r2 4 
+ (1 — 29)/ER, c— US — 29)/1241 —q)) 
(ahT, R/(k + hR)) r (92) 


varying from O at the center of the sphere to (93) at 
its surface 


u(R) = ((1 — 29)/E)t, R — ((3q — 4)/(12(1 — q)) 
(ahT, R2/(k + hR)) (93) 


From (52) - (87) - (88) and (38) - (39) - (40) the 
non-vanishing components of the Strain tensor are 


Sa tr) =((1— 29)/E)t, + (ahTo M(1201 =) (k + hR))) 
(2(1 + q)r — (9 — 2q)R) (94) 


Sp (r)=S (r)=((1— 19)/E)t, + (ahTo (Mk + 
= hRJCZO — q) (1 + qr — (5 — 29)R) (95) 


and using (41) - (42) - (43) the non-vanishing com- 
penents of the Srtess tensor are 


To (r)=t+ (aEhT, M12(k + hR) (1 q) (1 — 29))) 
((3 — g)r — (5 — 29)R) (96) 


Too(r) = Toolr) =t; + (Eh To /(k +hR)) ((5 — 
2q)/(12(1 — q))) (R +) (97) 


Therefore, the stress-state at the center of the 
sphere is pherical 


Ta (0) = Ta, (0) = Tt, (0) =t, + (15 — 29)/(12(1 
— q) (1 2q))) (ahT,R/(k + hR)) (98) 


and at the surface of the sphere is hemi-epherical 


TRI=t,+((8-— 3q)/(1201 — q) (1 29))) (aEhT, 
R/(k + hR)) (99) 


To(R)=To(R)=t;— U5— 29)/(611 — q) (1 
— 2q))) (aEhT, R/(k + hR) (100) 


Iv — HOLLOW SPHERE AND SPHERICAL LAYER 
IV.1 — Problem 4 


We consider an hollow sphere of inner and outer 
radiusec R, and R,, inner and outer surface pressure 
and temperature p,, T, and Pp, T respectivelly-The 
boundary conditions are 


T, (R,) — P; 
TiRg) = To 


TiR;)-p o THRj)=T, 
(101) - (102) - (103) - (104) 
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or, replacing (34) in (101) - (102), and (41) in (103) 
(104), and solving 


A=T,/R, —(R$/RyMT;R, — To R MRS — R$) 
(105) 


— pR2 p? ê — 83 3 
B=R2RZ(T, R$-ToR,DMR—R$—R)) 1106) 


C=(1-2q) p,/E- ((1 — 29)/E) (R$(p; — pj)/ 
NR$—R$)) + (a(3 — q)/(6(1 — q))). 


(RÊMR, —R5)MR$— R$ RMT, /R;—(T, R$— 
“To R,)/MR$— R$)) — (29R$/(0301 — q))) (107) 


D=((1+q)/2E) (R$R$(po — pi J/MR$— RT) 
US qt +aq)/(1201 — 0) (1 — 29))) (aRjR$ 
(Rj— R$) M(R$— R$)MT /R, — (R$/R, MT, Ro — 
— To Rj)/(R$ — R$)) + (al — q)/(301 — q)(1 — 29))) 
(aRÍRÍ(RZ— R2)/(R$ — RÍDUT, Ro — ToR|/R$- 
— R$)) (108) 


From (99) - (100) and (364), the temperature is 


T(r)=T,r/R, + R$(T, R$ ToR/ MR 
“RO ARZM TR) (109) 


and the non-vanishing component of heat flux is, from 
(99) - (100) - (36) 


Jr) =kT,/R, — (kR$ (TR, — T;R)) MR 
R$JJ(N/R, + 2R2/7º) (110) 


assuming the values at the surfaces of the sphere 


HR), =—kT;/R,—KRSUT,R9— To R,) MIR — 

— R$)MN/R,+2R$/R3) (111) 
MR$)=—KT4/R, —kR$UT;R$— To R, )MRZ— 
-“R$DMIR, + 2/R5) (112) 


We could obtain the displacement vector, and the 
Strain and Stress tensores replacing (105) - (106) 
(107) - (108) in (35), (38) - (39) - (40) and (41) - (42) 
(43); for example, the displacement vector is 


u(r) = (a(1— q)/(12(1 — q)MT,/R,— (R$/RjMT, 
R$ R$R,)MR$—R7)) 1 +[(1 — 29) p,/E — 

— ((1— 29)/E) (R$ (p; — py )/(R$ — R$)) +la(3 - 
-“q/(6(1 — q)MR$(R,— RIMR$ RIR) 
(T//R,— (R$/R,MT,R$— To R$) MRS — R$N) - 

— (ZagR$ (301 — q) (1 + Ro (R$ — R$JMRZ— R7)) 
UT,R, — T;R,)MR$ R$] r+ (alt +03 
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— qUMARI(T,R;— T;R, DMR R$) ro! + 

+ [((+ q)/2EMRÍRS(p,; — py) MRS — R$)) + ((3 — 
— qU/(12(L — q)(1 — 2q))) (aRJ R$ (R$ — R$)MRZ— 
-RÍIT,/R, —IRÍ RMT,R;—T;R,) MR$— 
— R$)) + (aq(1 + q)/(3(1 — q)(1 — 2a)))RÍR$ (R$ — 
— RIJMIRS — RÍDUT,Ro — To R,)/(R$— R$)) Jr? 


and the radial-radial componeut of the Stress tensor 


Tel) = (2E(3 — q)/(6(1 — q)(1 — 2q))) (T,/R, — 
—(R2/R,MT,R9— To RDMR$— RP r+ [py — 

— R$ (p,— PD/MR$— R$) + (E(3 — q)/(601 — q) — 
— 29))) (R$(R,— R9)MR$ RO — Rj) MT,/R, — 

— (R2/R MT, R$—To R,)/(R$ — R$)) — (2Eq/(301 — 
— q) (1 —2q))aRG(1 + R$(R$— R$)/(R$ — R$)) 
((T,R4— ToR,)/MR$— R$))] + (29E/(5(1 — 2q)))R$ 
R2(T,R,— T;R,DMRZ— R$) 2 [R$RS(p, — 

— po )/MR$— R$) + (E(3 — q)/(601 — q)(1 — 29))) 
(aRIRS(R,— R9)/MR$— R$MT,/R, —(R$/R,) 

(T;R9 — ToR,)/MR$— R$)) + (ZagE/(3(1 — q) (1 

— 2q))MRÍR$(RZ — R$JMRS — R7)) UT, Ro 

— To R,)MR$— R$] 1? 


IV.2 — PROBLEM 9 


We consider an hollow sphere of inner and outer ra- 
diuses R, and Ros pressuro and temperature p, and T, 
at the inner surface, and displacement Us and heat flux 
J, at the outer surface. Hence, the boundary con- 
ditions are 


TolR)=0, TiR,))=T, ulB,)=u, JMRg)=95 


(115) - (116) - (117) - (118) 
Reaplacing (41) e in (115), (34) in (116), (35) 
in (117) and (36) in (118), and solving for A, B, € 
and D, we obtain 
cida 3 3 : 
A=T//R,— (R$/kR,MUkT, + J;R,)/(R9 + 2R7)) 
B =(R$R$/kMKkT, + J,R,)/(R$ + 2R7) (2ERG/MT + 
+ q)R7) + E/((1 — 29)R7)) D (119) - (120) 
D=(p,R5— Eu,/( — 2q)) + (aER,/(6(1 — q)))MU3 — 
— q)R,/1 — 29) — (1 + q)R9/2MT,/R, — (R$/kR ,) 
(KT, + J9R,)/(R$ + 2R7)) + (aE/(3(1 — q(1 — 
— 2q)))29R5 /R$— (1 +a)/R9 JUR, R$/kMKT, + 
Jo RjJ/(R$ + 2R7)) (121) 


C=(1-2a)p1/E-— (aR,(3-q)/((11 — QT, 
[Rj — (R$/KR MKT, + Jo R$) MRS +2R$)) — 
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— (2q/(3(1 — q))) (aR, R$/KMMKT, + Jo R$ )/MRS + 

+ 2R7)) + (2(1 — 2q)/(1 + q))(D/R7) (122) 
Replacing (119)-(120)-(121)-(122) in equa- 

tions (34) to (43), the temperature, displacement, heat 


flux, strain and stress tensores, colud be obtained, as 
in the proceding problems. 


IV.3 — PROBLEM 6 


We consider an hollow sphere of inner and outer 
radiuses R | and R5 , displecement u, ant and heat J| 
through the inner surface, tangential 'stress t, and heat 
convection from a fluid at free temperature pé at the 
outer surface, The boundary conditions are 


ulR)=u,  JHRj)=3) TIR)=TiRd=t 
JHR9)=h(TIR;3-To) 


(123) - (124) - (125) - (126) 


Replacing (35) in (123), (36) in (124), (42) in 
(125) and (24) - (36) in (126), and solving for A, B, 
C, and D, we obtain 


fi 3 
A=-—J/k+2R$(khT, + J, (hR$+k))/ 


M2kKRÉ(hRo + k) — R$ (2k — hR$)) (127) 


B=RSR$(khTo + 9, (hRo + K)/ZKRZ (hR$ +) - 
— kR$(2Zk — hR$)) (128) 


C=u,/R, — (a(1 + q)R$/(1201 — q) (— J,/k+ 

= 2R$ (khT, + dj ((hR5 + k)) MZKRÊ (hR5 + k) — kR$ 
(2k— hR$ DD) — (1 +31 — q))) Rj R$ (khT, 

+ 3, (hR$ + k)/(2kKR$ (hR$ + k) — kR$ (2k — 


— hR$))) — D/Fj (129) 


D=((1+aq)(1 — 29)/3q)(R7/E) (Eu, /(1 — 29) — ty 
Ry) — (aRj(1 + q)(4 — 3q)/(369(1 — q)))( — 3, /k + 
= 2R3 (khT + J, (R$ + k)) /(ZKR$ (k + hR$) — kR$ 
(2k — hR$))) — aRf R$ ((1 + q)(3a — 1)/(9a(1 — q))) 
(khTo + Jj (R$ + K))/(2kR$ (hR 9 + k) — kR9 


(2k — hR$))) (130) 


Replacing (127) - (128) - (129) - (130) in eque- 
tiois (34) to (43) we ?ould obtain the temperature, dis- 
placement, heat flux, strain and stress tensores., 


IV.4 — PROBLEM 7 


We consider a spherical layer of radius R and 
thickness e, at inner and outer surface temperature and 
pressure p,, T, and ps, Ta. The boundary condi- 
tions are 
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T(R— 6/2) =T, T(R+e/2)=T; TU(R—- 0/2) =p, 
T,(R+e/2) =p, (131) - (132) - (133) - (134) 


Replacing (34) in (131)-(132) and (41) in (133) 
(132) and (41) in (133) - (134), and solving for A, B, 
C and D, or, directely from equations (105) - (106) 
(107) - (108) of problem 4, using the relations (40), 
— that make boundary conditions identical — , and ne- 
glecting higher powers of e, we obtain 


B=((T|—7T;)/30) Rº 
(135) - (136) 


A=T,/R+(T —T,)/3 


C=(1-2q)p,/E — ((1 — 2q)/E)(Ripo — p, 1/30) + 
(2aR(3 — q)/(18(1 — q))) (T,/R+(Ta —T, )/36) 
— (10agR/(9(1 — q))JM(T, — T5)/6) (137) 


D=((1+q)/2E)Rº (po — p,)/30) + (aRº/3)((3 — 
qU1 + q/(1201 — q) (1 — 29) )MT,/R — (T, — 
—T,)/36) + (aq(1 + q)/(5(1 — q)(1 — 29) (2Rº (T, — 

— To )/9) (138) 


From (125) - (126) and (34), the temperature is 


Tt) =T,/R+MT,—T5)/30)lr — R3/r? (139) 


and from (125) - 136) and (36) the non-venishing com- 
ponent of the heat flux is 


Jtr)=kT,/R+K(T, —T5)/30 + 2k(MT, — To)/ 


/36)(Rº/rº) (140) 


We could obtain the displacement vector, strain 
and stress tensores, replacing equations (135) to (138) 
in equations (34) and (38) to (43). For example, the 
displacement vector is 


ulr) = (a(1 +aq)/(1201 — q)))(T,/R-—T, 

—T9)/360)r2 + [(1 — 29)p, /E — ((1 — 293/E) 

(R(pz — p/)/363 + (2a(3 — q)/(18(1 — q RIT,/ 

/[R+44To—T,)/30) — 10a9R/(3(1 — q))) UT, — 

— T2)/30)] r + (alt + a) (301 — q) RÉT, — 

— T9)/38) r”! + [((1 + q)/2E)R* (lp, — p,)/30) + 

- (a(3 — qu(1 + q)/(36(1 — q) (1 — 29) DR“ (T,/R — 
(Ty — T, 1/36) + (2aq(1 + q)/(9(1 — q)(1 — 29))) 

Rº((T, — T,)/30)] ro? (141) 

and the radial-radial component of the stress tensor 

Ta (1) = (aE(3 — q)/0(1 — qu(1 — 29))M(T,/R — 

-(T,—T,)/30) r+ [py — Rip, —p,)/30 — (E(3 — 

— q/(6(1 — q)(1 — 2q)))(2e/3)M(T,/R— (T, — T,)/ 


[38) — (43Eg/(9(1 — q)(1 — 29)) (UT, — T,)/30) | + 
+ (2aE/(3(1 — 2q)))Rº ((T) — T4)/36) 172 — 
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—[Rº(p, —p,)/3e + (aE(3 — q)/((18(1 — Q) (1- 
— 2q)))R“(T)/R — (T|— T,)/36) + (4Eag/(9(1 — 29))) 
ReMTy=—T,M30)] (142) 


V — SPHERICAL CAVITY IN INFINITE MEDIUM 
V.1 — Problem 8 


We consider a spherical cavity of radius R, in an 
infinite medium, the pressure and temperature at the 
surface of the cavity being p and T, respectivell. 

We remark that, in all spherical cavity problems, 
since the temperature (34) and Displacement (36) must 
remain finite as r tends to infinity 


(143) 


The boundary conditions of this problem are 


TR) =To, T;(R)=p (144) - (145) 


Replacing (34) in (144) and in (145), using (143), 
and solving for B and D 


B=T,R?D=-((1+q)/2E)pRº+ (a(1 + q)/(3(1 — 
— qu1 — 29) T;Rº (146) - (147) 


From (143) - (146), and (34), (36)), the tempe- 
rature and non-vanishing component of heat flux are 


J(r)=2kT,R2/rº J(R)=2kT,/R 
(148) - (149) - (150) 


Ttr) =T5R2/r? 


the former decreasing steadily from its value (144) 
at the cavity's surface, and also the latter, whose value 
at the cavity's surface is given by (150). 

From (143) - (146) - (147) and (35) the displa- 
cement vector is 
u(r) =— (( + q)/E)pr”/r? + (aT5(1 + q) /(301 — 


qUMR2Z/r + (q/(1 + q))R3/r2) (151) 


so that the displacement of the cavity's surface is 


u(R) = — ((1 + q)/2E)pR + aTÇR(1 + 29)/(3(1— 29)) 
(152) 


From (143) - (146) - (147) and (38) - (39) - (40), 
the non-vanishing components of the strain tensor are 


Sr) = (1 + q)/E)pRê/ré — (97, (1 + q)/ (301 — 


— q))M(R2/r? + (29/(1 — q))R“/rº) (153) 
Soy (1) = Soap lr) = — (11 + q)/2E)pRÉ/r+-. (aT5 (1 + 
+ q)/(3(1 — q))MR2/r? + (q/(1 — 2q))R3/rº (154) 
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and from (41) - (42) - (43) the non-vanishing com- 
ponents of the stress tensor are 


To (1) = pRÊ/rº + (29EaT, (1 — q) (1 — 2q)) DR? / 


fr? — R$/rº) (155) 
Toolt) = Top (rn) = — (p/2)Rº/rê + (28ET5/(5(1 — 
— 29) MR 2/12 = (1/(1 — q))R$/rº) (156) 


The stress-state at the surface of the cavity is 

hemi-spherical, specified by (145) and 
To (R) = Ty (R) = — p/2 — 29EaT /(3M1 — 
— q)(1 — 2q)) (157) 


V.2 — Problem 9 


We consider a spherical cavity of radius R in an 
infinite medium with displacement u, and heat flux 


Jo through its surface. The boundary conditions are 


u(R)=u, JMHR)=J, (158) - (159) 


Replacing (35) in and (36) in (159)( using (143), 
and solving for B and D 


B=J,Rº/2k D=u,R?2-((L=aq)/(3(1— 
— q)JaJ,R*/2k (160) - (161) 


From (143) -(160) and (34) the temperature is 
given by (162), and decreases steadily from its value 
(163) at the surface of the cavity. The non-vanishing 
component of the heat-flux is, from (145) - (160) and 
(36), given by (164), and also decreases steadily from 
its value (159) at the surface of the cavity, 


T(r) = (J,/2k)R3/r? T(R) =J,R/2k Jr) = 3, R%/rº 
(162) - (163) - (164) 


From (143) - (160) - (161) and (35) the displa- 
cemente vector is given 


u(r) =u,R?/rº + (a(1 + q)/(301 —a) (J,R?/ 


/[2kMR/r — R?/r*) (165) 


its value at the surface of the cavity is (158), and its 
maximum is at tha radius Ra 


Rm = 2R-(3(1— q)/(1 + q))(2k/ad, R) 


Rá UR) =u,R2-RÍ+((L+aq)/(3(1 — q))) 


(aJ, R$/2k) (167) 


From (143) - (160) - (161) and (38) - (39) - (40), 
the non-vanishing components of the strain tensor are 


Se ll=-—2u,R2/r2+((1 +q)/(3(1 — q))(as,R/ 
/2k) (2R%/7º — R9/r2) (168) 
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Sol) =Sopl)=u,R2/02+ 1 + q)/(31 — 


— q)))aJ 5, R/2kMR2/r2 — R$/r3) (169) 


and from (41) - (42) - (43) the non-vanishing of the 
stress tensor are 


To tr) = — (2E/(1 + qlluç R2/rº + (2E/(3(1 — 


— q)))aJ, R/2k)(q/(1 — q))R 2/12 — R2/1º) (170) 


Too lr) = Toy (r) = (E/(1 + q)JuçR2/rê + (E/(3M1 — 
— q))MaJ,R/2k)R?/r3 + (2E/(3(1 — 2q))) (ad, R/ 
/2k)R Z/r* (171) 


The stress-state at the surface of the cavity is 
hemi-spherical 


T, tr) = — (2E/(1 + qlju, /R + (2E(3q — 1/30 — 
— qh(1 — 2q)))aJ, R/Zk) (172) 

To (1) =To (r) = (EM + q)luç/R + (E(3 — 

— 49) /(3(1 — alli — 29)))aJ,R/2k) (173) 


V.3 — Problem 10 


We consider a spherical cavity of radius R in an 
infinite medium, with tangential stress t, at its surface, 
and heat convection from a fluid at free temperature T ã 
through its surface. The boundary conditions are 


T()=TI(R)=T, J(R) = h(T(R) — To) 


(174) - (175) 


Replacing (42) in (174), and (34) - (36) in (175), 
using (-43), and solving for B and D, we obtain 


=hT,RÍ/2k+hR) D=((1+q)/ERçRº— 
— (2(1 + q)/(311 — 29)))lahT, RÉ/(2k + hR)) (176) 


From (176) - (143) and (34), the temperature is 
given by (178), and (177) decreases steadily from its 
value (179) at the surface of the cavity 


T(r) = (hT,/12k + hR$)R/rS 


T(R) = hT5 R/(2k + hR) (178) - (179) 


From (176), (143) and (36), the non-vanishing 
component of the heat flux is (180), and it decreases 
steadily from the value (181) at the surface of the 
cavity 


J(r) = (2khT,/(2a + hR)JR$/rº 


J(R) = 2khT, /(2k + hR) (180) - (181) 


From (143), (176), (177) and (35) the displerce- 
ment vector is 
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ulr) = ((1 + q)/(3(1 — q)))(ahTo R/(2k + hR)) Rº/rA 
+ (UL — q)/Edt, RÊ/r? — (201 + q) 31 — 
— 2q)))ahT, R/(2k + hR)JR?/r? (182) 


and the displacement of the cavity's surface is 


u(R)=((1+q)/E,R—((1 + q) — q) — 
— 2q)))ahT,, Rº/(2k + hR)) (183) 


From (143) - (176) - (177) and (38) - (39) - (40), 
the non-vanishing components of the Strain Tensor are 


Se ()=—(U1+a)/(301 —q))) (ahTo R/(2k + 
= hRDJR2/r? — (2(1 + q)/EXÇR$/ ré + (4(1 4 
+ q)/(301 — 2q)))ahT, R/(2k + hR)JRÊ/rº (184) 


S (M)=S ()=(1+q)/(301 — q)))ahT, R/(Zk + 
E hRDR2/r2 + (1 + q)/EN,RÊVró + (21 4 
» q)/43(1 — 2q))) (ahT, R/(2k + hR)JRÊ/rê (185) 


and from (41) - (42) - (43) the non-vanishing compo- 
nents of the Stress Tensor are 


Ta (rn) = — 2t,R?/0º+ (2E/(3(1 — 29))))ahTR/(2k + 
+ AR) (ZRZ/rº + (q/(1 — q))R?/ 12) (186) 


Til=T()=t,Rê/rº + (2E/(3(1 — 2q)))ahT R/ 
2k + hR)) (R$/r2 - R3/73) (187) 


The Stress-state at the cavity surface, is hemi- 
spherical, specified by (174) and 


T (R)=—12, + (2E(2 — q)/(301 — q) 
(1 — q)))(abT, RJ/(2k + hR)) (188) 
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FUNÇÃO APROVISIONAMENTO — GESTÃO DE STOCKS 
SUA INFLUÊNCIA NA ECONOMIA DAS EMPRESAS E DOS PAÍSES 


RESUMO 


Justifica-se teoricamente das vantagens económicas 
da existência dum serviço único de Aprovisionamentos 
dentro duma empreza, ou do agrupamento de empre- 
zas para um Serviço de Aprovisionamento comum, ou 
ainda, ao nível comercial dum país, as vantagens do 
recurso ao comerciante/representante/distribuidor em 
vez de importarem. 


INTRODUÇÃO 


Recentes disposições legislativas (1) disciplinaram, 
em Angola, o regime das importações. Tal legislação 
implica ainda que no nosso país se passa a tomar em 
atenção a importância que as técnicas de APROVISIONA- 
MENTO, nomeadamente a Gestão de Stoks, têm na 
economia das empresas e dos países. 

Há já alguns anos que, por força da minha vida 
profissional, tive de inteirar-me do que era a função 
Aprovisionamento e, actualmente, encontro-me mesmo a 
organizar tal serviço numa importante empresa Ango- 
lana. Daí o sentir com maior acuidade a importância de 
tal legislação e o que ela pode representar para o 
desenvolvimento do país. 

Transcrevemos, da Circular n.º 22 da Direcção dos 
Serviços de Comércio da Província de Angola, o terceiro 
parágrafo da página 7: 


— «Sempre que possível, devem os importadores 
inscritos na classe única recorrer ao Comércio local 
para aquisição do que necessitem para as suas activi- 
dades, evitando o recurso directo à importação, que é 
muitas vezes anti-económica, o que contribuirá para o 
aparecimento de comércio especializado de apoio à 
indústria.» 


Exprime tal afirmação o desejo da governação para 
que os importadores inscritos na classe única limitem 
a sua actividade de Importação comprando em Angola, 
aos Representantes dos Fabricantes, os produtos de que 
necessitam, 


FRANCISCO PRISTA CAETANO 


SYNOPSIS 


This is a theoretical justification of the economic 
advantages which results when, within one company all 
buyng and stocks are concentrated in only one service 
and, simmilarly, the same happens when some Compa- 
nies came to an agreement for a unic Procurement Ser- 
vice, or at the country commercial level, the companies 
decide to buy from the dealers instead of importing 


Este procedimento permitirá dentro de algum tem- 
po, aos Organismos Económicos, ter uma ideia real dos 
consumos de Angola e, assim, orientar a criação de 
parque Industrial com bases realísticas. 

Além desta acção de futuro há ainda uma resul- 
tante económica imediata que será tanto mais especta- 
cular e efectiva quanto mais rapidamente os senhores 
Industriais-Importadores resolverem comprar aos repre- 
sentantes em vez de importar, e exigir destes a consti- 
tuição de stocks capazes de satisfazer as necessidades 
das suas indústrias. 

Foi o ver «oficializado», um ponto de vista que há 
tanto tempo tenho vindo a advogar, que me levou a 
fazer algumas considerações teóricas simples mas que, 
no entanto, provam a justeza da legislação que citei como 
tendente a uma economia imediata, não só para o pró- 
prio país como também para as empresas. 

O assunto Gestão de Stocks tem sido tratado por 
vários estudiosos sendo Wilson um dos mais conhecidos, 
e é a partir dos seus estudos que faremos as nossas 
considerações. 


GESTÃO DE STOCKS 


A quantidade económica Q a encomendar de cada 
vez é dada pela fórmula de Wilson. 


“4 


na qual é 


(1) Decreto-Lei N.º 478/71, de 1971. Nov. 06, Portaria Ministerial N.º 703/71, de 1971. Dezembro. 15 e Decreto N.º 552/71 


da mesma data. 
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A — Custo unitário de efectivação duma encomenda, 
N — Consumo anual em unidades. 
P — Preço unitário pago ao fornecedor. 
r — Coeficiente de custo de posse de stock. 
(0.15 a 0.20 do valor do stock). 


A fórmula de Wilson obtém-se a partir da consi- 
deração do custo total C de Aprovisionamento, o qual 
é constituído por: 


— Custo de aquisição. C1=N.p 

N 
— Custo da encomenda. C2=A- E 

Q r.Ã 
-— Custo de posse de stock C3 =r.p.. e + E 


C=C1+4C2+C3- 


N Q r.Ã 
=Np+rA —— + tb. - E 
a a ii 


em que N, p, À, O e r são as mesmas entidades indi- 
cadas atrás. 

Consideremos n empresas que teriam de Aprovi- 
sionar um mesmo bem de consumo e que o podem 
fazer independentemente ou agruparem-se para um 
Aprovisionamento único. 

A quantidade total a encomendar de cada vez pela 
empresa seria 


n 
Q, - *X Qn=n Om 
l 
n 
em que Om = 1 Qn seria a quantidade média a en- 


n 


comendar por cada empresa, 
Do mesmo modo, o consumo anual total N, das n 
empresas, seria 


e 
em que N, — 71 Nn é o consumo médio anual de 


n 


cada uma das n empresas. 
Portanto, se uma empresa única comprar para as n 
empresas, o custo total será: 


Nt Q r.Ã 
CN. PFA— + rp + — 
ati fds a, "a 2 


em que O, será a quantidade económica a comprar de 
cada vez pela empresa única. 
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Pela mesma dedução de Wilson será: 


3 2xAxN, 
= 
ou 
[2xAxnxN C2xAxN 
X Xn x a / X X 
Ei Vs L =y a «V——— 
F.p Fr.p 
mas 
/2xAxN a: 
W= — = Om = 
r.p n 


que representa a quantidade económica média que teria 
de ser adquirida por cada uma das empresas, indivi- 
dualmente. 

Podemos portanto concluir que 


ou 


OG -— x0,= x 


pa t — 


V n VY n 


mo 4 53 
|) 


Isto é, a quantidade económica a adquirir por uma 
só empresa é uma fracção da totalidade das quantida- 
des económicas a adquirir por cada uma das empresas 
e essa fracção é exactamente o inverso da raiz quadrada 
do número de empresas que se agruparem. 

Apresentamos em seguida um exemplo valorizado 
com que desejamos ilustrar o que dissemos. 

Vamos considerar um modelo constituído por n = 5 
empresas industriais. 

Consideremos um bem de consumo necessário e 
fundamental à vida daquelas empresas. 

Para facilidade de exposição consideraremos que 
todas aquelas empresas teriam o mesmo valor de con- 
sumo anual e que todas têm bem organizados os seus 
serviços de Aprovisionamentos, com os respectivos sec- 
tores de Armazém, Compras, Gestão de Stock, Recepção 
e Técnica e que, portanto, têm os mesmos valores carac- 
terísticos para: 


= 2000800 (encomenda ao estrangeiro). 
10 000/ano. 

= 20800. 

= 0.2 (20%). 


Â 
N 
p 
r 


/2 x2 000 x 10 000 
o Gado “AH 


/ 2AN. 


o=V “5 =4 


A quantidade total que todas as 5 empresas teriam 
de comprar seria portanto: 


O, =5 x Q = 15800 


t 
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A quantidade económica que a empresa única com- 
praria seria: 


x 0 =. |! s 15800 =7085 


vY 6 


Se considerarmos os valores 


OQ, = 15 800 x 20800 = 316 000800 
O, = 7085 x 20500 = 141 700$00 
Portanto: 


Economia somente em empate de capital174 300800. 


É interessante notar que este raciocínio também se 
aplica aos vários sectores utilizadores dentro duma 
mesma empresa e, portanto, podemos concluir e deve- 
mos recomendar que não devem os vários sectores du- 
ma mesma empresa constituir «stocks» particulares mas 
sim requisitarem ao Serviço de Aprovisionamento 
somente os materiais de que necessitem para uso 
imediato. 

O Serviço de Aprovisionamento comportar-se-á 
como a empresa única a comprar e «stockar» para todos 
os sectores consumidores da empresa. 
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“MUNDO TÉCNICO” 


MATERIAL PARA SEGURANÇA INDUSTRIAL 


Foi lançado no mercado um material transparente que pode constituir uma protecção muito eficaz para O 
pessoal de oficinas de soldadura, mecanização, moldagem, carpintaria e fundição. 

Este material é tão resistente que pode impedir a passagem de uma bala de 5 mm disparada a 23 m de 
distância. 

Trata-se de uma lâmina de acetato de celulose, que é inestilhaçável e adequada para a fabricação de 
protecções de máquinas, óculos de protecção, visores e janelas de segurança. 

A sua resistência ao impacto é, pelo menos, 12 vezes maior que a chapa acrílica e duas vezes e meia 
mais rígida que a de cloreto de polivinilo. 

Ainda que seja perfurado o material por fragmentos a grande velocidade não se estilhaça. 

O material reune boas propriedades priorresistentes e um ritmo de combustão que foi classificado de lento 
a autoextintor. 

Este acetato de celulose é de fácil moldagem, podendo ser moldado normalmente. 

Não se amarela com o tempo, e vende-se em várias cores para melhorar a protecção da vista na soldadura. 
(DYNNA) 

L. A. Petri. 

S. A. Gerona 78 Barcelona-9. 


RENOVAÇÃO DO REVESTIMENTO DE UM FORNO DE REVENIR COM FIBRA 
CERÂMICA 


Num forno de revenir da «Safir», França, tirou-se recentemente o revestimento inicial de ladrilhos refrac- 
tários e colocou-se no seu lugar outro de fibra cerâmica, em 70 horas, depois das quais o forno se achava em 
condições de começar a funcionar imediatamente. 

À renovação do revestimento com materiais refractários normais teria exigido muito mais tempo e o período 
de intercepção na produção teria de ser prolongado enquanto não terminasse de secar. 

A fibra cerâmica utilizada foi Triton Kaowool, material apropriado para emprego até 1260º C e que pesa a sexta 
parte do ladrilho refractário isolante. 

É completamente elástico e insensível aos efeitos da expansão técnica, porque a expansão de cada uma 
das fibras fica absorvida dentro do material. 

No sitio mais denso deste revestimento a condutividade térmica é inferior a 35% da do ladrilho refractário 
isolante e a sua capacidade de acumulação de calor é '/, aproximadamente desse ladrilho. 

Estas propriedades permitem obter uma redução considerável nos tempos de caldeação e de esfriamento 
e a economia no consumo, chega a ser de 15%. (DYNNA) 


BOMBA PNEUMATICA 


A ENERPAC apresenta uma nova linha de bombas pneumáticas de alta pressão denominadas «Bomba pneu- 
mática dupla». 

Esta bomba possui as seguintes dimensões 380xX320x308 mm, completamente portátil, pode ser ligada a 
qualquer linha de ar comprimido da oficina, incluindo o compressor móvel. 

A máxima pressão hidráulica é de 700 bars (ou 10 000 Psi) e pode ser regulada pelo ajustador da linha de ar. 

A originalidade desta bomba pneumática reside na aplicação de 2 pequenos motores pneumáticos: 

— Um para o avanço rápido, admitindo o alto fluxo e baixa pressão da linha de ar a 7 kg/cm2 (1 galão 
por minuto com 100 Psi — pressão pneumática). 

— Um para alta pressão, permitindo pequeno fluxo de óleo a alta tensão, 

(0.16 I/m a 700 kg/cm2 (pressão pneumática 7 kg/cm2) 10 in3 a 10000 Psi (pressão-pneumática 
100 Psi). 

Esta bomba tem um ruído inferior a 80 decibeis (ROMAR). 
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ARMANDO NUNES PIRES CAMEIRA 
ALGUNS CONCEITOS BÁSICOS 
RELATIVOS A PROBLEMATICA DOS TRANSPORTES 


Técnica N.º 419 — XLVII — 1-1973, p. 159 - 164 


C. D. U. 656.01 


O presente artigo contém alguns conceitos gerais sobre a 
engenharia de transportes. Constitui o texto de uma lição 
dada num curso de férias aos alunos de eng.” civil, promo- 
vido pelo Departamento de Comunicações do |. S. T. 


OLIVÉRIO SOARES C. D. U. 537.311.6:621.3.126 


CONSIDERAÇÕES SOBRE A INDUTÂNCIA 
DUMA MAQUINA DE CORRENTE CONTINUA 


Técnica N.º 419 — XLVIII — 1-1973, p. 165-170 


O autor expõe o conceito de indutância, e aplica-o às má- 
quinas de corrente contínua. 

O artigo é a compilação de ideias, para desenvolvimento 
ulterior, em que são invocadas relações importantes. 


C. D. U. 621.396:551.596.7 
JOSE A. C. SARAIVA MENDES 


MEDIDAS DE ATENUAÇÃO A 11 GHz 
EM AMBIENTE CHUVOSO 


Técnica N.º 419 — XLVIII — 1-1973, p. 171-179 


Este trabalho destina-se a apresentar alguns resultados 
experimentais provenientes de medidas de atenuação, feitas 
em dois percursos a 11 GHz durante os anos de 1970 e 
1971. Analisa-se também um método de cálculo proposto 
por P, Quarta, permitindo o cálculo da distribuição estatis- 
tica da atenuação instantânea prevista num percurso, a 
partir da distribuição estatística da precipitação média 
horária relativa a qualquer ponto da região onde se encontra 
esse percurso. 


D. F. G. DURÃO 
INVESTIGAÇÕES EXPERIMENTAIS EM JACTOS COAXIAIS 
Técnica N.º 419 — XLVIII — 1-1973, p. 181-189 


Apresentam-se medidas da média e flutuações de veloci- 
dade assim como das flutuações do sinal eléctrico do fio 
quente. Elas são uma investigação preliminar do escoa- 
mento para confirmar a possibilidade de ser usado em me- 
dições detalhadas dos parametros da turbulência. A concen- 
tricidade do escoamento é confirmada e tem um perfil de 
velocidade média totalmente desenvolvido a 30 diâmetros. 


JAIME M. COSTA OLIVEIRA 
CARLOS A. MORAIS DA VEIGA 


DETERMINAÇÃO EXPERIMENTAL DA FUNÇÃO 
DE TRANSFERÊNCIA DO REACTOR PORTUGUES 
DE INVESTIGAÇÃO A POTÊNCIA ZERO 


Técnica N.º 419 — XLVIIH — 1-1973, p. 191-200 


Mediu-se a função de transferência do Reactor Português 
de Investigação (RPI), a baixa potência, pelo método de 
modulação da reactividade, na banda de frequências com- 
preendida entre 5.10 ——? Hz e 65 Hz. Descrevem-se as 
caracteristicas dos dispositivos experimentais e da técnica 
de medição utilizados, bem como o princípio do tratamento 
numérico dos dados experimentais. 

OQ modulador concebido para esta experiência possuía boas 
características de estabilidade para as frequências estudadas 
o a técnica adoptada na detecção de análise da resposta 
do reactor revelou-se muito eficaz. Os resultados obtidos 
são comparados com os fornecidos pela teoria cinética 
pontual, constatando-se a influência dos harmónicos espa- 
ciais mo comportamento dinâmico do reactor, a qual não 
é tida em conta pelo modelo pontual, Dos valores da função 
de transferência deduz-se para o tempo de geração dos 
neutrões instantâneos do RPl o valor , = (56+2).10-ºs 


para uma fracção efectiva dos neutrões retardados Be 
= (760 + 23). 10 —+. 


C. D. U, 621.039,5 


LUIS MANUEL BRAGA DA COSTA CAMPOS 


ALGUNS PROBLEMAS DE SIMETRIA ESFÉRICA 
NO CONTORNO (FRONTEIRA) EM REGIME PERMANENTE 
NUM MEIO ANISOTÉRMICO E ELASTICO 


Técnica N.º 419 — XLVIIH — 1-1973, p. 201-210 


O objectivo deste artigo é estudar uma classe de problemas 
inter-relacionados, respeitantes à determinação da tempe- 
ratura, vector de fluxo de calor, vector de deslocamento, 
tensor de deformação, tensor de tensão e vector de tensão, 
de um meio linear, homogéneo, isótropo e elástico, em 
regime permanente, com condução de calor, para 4 geome- 
trias esféricas — esfera oca, camada esférica, esfera maciça 
& cavidade esférica num meio infinito —, e condições na 
fronteira especificando a deslocação superficial, a pressão 
ou a tensão tanjencial, e a temperatura, o fluxo de calor 
ou condições de convecção. 

Alguns problemas de simetria esférica às condições na 
ed em regime permanente num meio anisotérmico 
o elástico. 


FRANCISCO PRISTA CAETANO C. D. U. 658,7 


FUNÇÃO APROVISIONAMENTO —- GESTAO DE STOCKS 
SUA INFLUÊNCIA NA ECONOMIA 
DAS EMPRESAS E DOS PAISES 


Técnica N.º 419 — XLVIIH — 1-1973, p. 211-213 


Justifica-se teoricamente das vantagens económicas da 
existência dum serviço único de aprovisionamentos dentro 
duma empresa, ou do agrupamento de empresas para um 
Serviço de Aprovisionamento comum, ou ainda, ao nível 
comercial dum país, as vantagens do recurso ao comer- 
ciante/representante/distribuidor em vez de importarem. 
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Synopsis of articles published in «Técnica» n.º 419 
XLVIH — January 1973 


JAIME M. COSTA OLIVEIRA 
CARLOS A. MORAIS DA VEIGA 


ZERO POWER TRANSFER FUNCTION MEASUREMENT 
OF THE PORTUGUESE RESEARCHE REACTOR (RPI) 


Técnica N.º 419 — XLVII — 1-1973, p. 191-200 


U. D. C. 621.039.5 


The zero power transfer function of the Portuguese Research 
Reactor (RPI) was measured using the reactivity oscillation 
technique cver the frequency range of 5.10 --3 Hz to 65 
Hz. The experimental set-up and measurement technique 
used are described, as well as the numerical handling of 
the data. 

The reactivity oscillator designed for the present experi- 
ment had good stability at the frequencies studied and the 
technique used in the detection and analysis of the reactor 
output proved to be very efficient. The results obtained are 
compared to the point-reactor kinetics model and the effect 
of the space harmonics on reactor dynamics is observed. 
From the measured transfer function a prompt neutron 
generation time of A =(56+ 2). 10 -º sec is obtained 


with an effective delayed neutron fraction Sp = (760 + 23). 
10 —3. 


U. D. C. 539.3:536.24 
LUIS MANUEL BRAGA DA COSTA CAMPOS 


SOME SPHERICAL-SIMETRY BOUNDARY-VALUE 
PROBLEMS IN STEADY-STATE ANISOTHERMAL 
ELASTIC MEDIUM 


Técnica N.º 419 — XLVII — 1-1973, p. 201-210 


The purpose of this article is to study a class of interre- 
lated problems, concerning the determination of the tempe- 
rature, heat flux vector, displacement, vector, strain tensor 
stress tensor, and stress vector, of a linear, homogeneous, 
isotropic elastic msdium, in steady-state, with heat conduc- 
tion, for 4 spherical geometries — hollow sphere, spherical 
layer, solid sphere and spherical cavity In infinite medium—, 
and a range of boundary conditions, — given surface displa- 
cement, pressure or tangential stress, and temperature, heat 
flux, or convection conditions. 


FRANCISCO PRISTA CAETANO U. D, C. 658.7 


IMPORTANCE OF THE FUNCTION SUPLY 
— MANAGEMENT OF STOCKS 


Técnica N.º 419 — XLVIII — 1-1973, p. 211-213 


This is a theoretical justification of the economic advan- 
tages which results when, within one company all buyng 
and stocks are concentrated in only one service and, simmi- 
larly, the same happens when some Companies came to 
an agreement for a unic Procurement Service, or at the 
country commercial level, the companies decide to buy 
from the dealers instaad of importing. 
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ARMANDO NUNES PIRES CAMEIRA U. D. C. 656.01 


SEVERAL BASIC CONCEPTS ON TRANSPORTATION 
ENGINEERING 


Técnica N.º 419 — XLVIII — 1-1973, p. 159-164 


This paper contains several and general concepts on trans- 
portation engineering. It was a lesson given during holyday's 
course at engineer's students, promoted by Departamento 
de Comunicações of |. S. T. 


OLIVÉRIO SOARES U. D,. C, 537.311.6:621.3.126 


SOME CONSIDERATIONS ON THE INDUTANCE 
OF A DIRECT CURRENT MACHINE 


Técnica N.º 419 — XLVIII — 1-1973, p. 165-170 


The author explains the conceit of inductance and applies 
it to direct current machines. 

The article is the compilation of ideas for an ulterior 
development, in which important relations are invoked. 


U. D. C, 621.396:551.596.7 
JOSE A. C. SARAIVA MENDES 


ATTENUATION MEASUREMENTS AT 11 GHz 
UNDER RAIN CONDITIONS 


Técnica N.º 419 — XLVII — 1-1973, p. 171-179 


This paper shows some experimental results obtained 
during 1970 and 1971 on two paths at 11 GHz. The ana- 
lysis of a method proposed by P. Quarta is also carried 
out. Such method allows the evaluation of the statistical 
distribution of the instantaneous attenuation on a path, from 
the knowledge of the statistical distribution of the hourly 
mean rain rate measured at any point of the country 
where the path in located. 


D. F. G. DURÃO U. D. C. 532.544 
EXPERIMENTAL INVESTIGATIONS OF CO-AXIAL JETS 
Técnica N.º 419 — XLVIINI — 1-1972, p. 181-189 


Measurements of mean and rms velocity and rms electrical 
signal from a hot-wire are reported, These measurements 
represent the preliminary flow investigation necessary to 
confirm that the flow is suitable for detailed turbulence 
measurements. The concentricity of the flow is confirmed 
and the flows, as a result of the fully developed initial 
profile, show to have a fully developed mean velocity profile 
at 30 diameters from the jet exit. 
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Empresa de Sondagens e Fundações 
TEIXEIRA DUARTE, LL.” 


JBLICA, 42, 8º 
“Tu. Iodo é is LISBOA 


SONDAGENS GEOLÓGICAS 
CAPTAÇÃO DE AGUAS SUBTERRÂNEAS 
CONSOLIDAÇÃO DO SOLO 
INJECÇÕES DE CIMENTO 


FUNDAÇÕES DE TODOS OS TIPOS 
CAYES E TUNEIS 
BARRAGENS E PORTOS 


(UM QUARTO DE SÉCULO DE ESPECIALIZAÇÃO TÉCNICA) 


ARRANQUE INSTANTÂNEO 
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COM BATERIAS - PILHAS SECAS 


TUDOR 


novas instalações de vendas e assistência técnica: 
em LISBOA- R.Alferes Malheiro,5A-B Tel.7306 81 (A averasiL) 
no PORTO- Av. da Boavista, 7146-754 Tel. 694049 


